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SUMMARY 


A transfer function method, for predicting the dynamic 
responses of gear systems with more than one gear mesh is de- 
veloped and applied to the NASA Lewis four-square gear fatigue 
test apparatus. Methods for computing bearing-support force 
spectra and temporal histories of the total force transmitted 
by a gear mesh, the force transmitted by a single pair of 
teeth, and the maximum root stress in a single tooth are de- 
veloped. Dynamic effects arising from other gear meshes in 
the system are included. 

It is shown that the above response metrics depend funda- 
mentally on the dynamic transmission error of the mesh under 
consideration. The Fourier transform of the dynamic transmis- 
sion error of a mesh is expressed as a linear combination of 
the Fourier transforms of the static transmission errors of 
all meshes in the system. The complex frequency-dependent co- 
efficients in this expression, which we call transmission error 
dynamic increment functions, characterize the dynamic properties 
of the gearing system. Expressions for these functions are de- 
rived and plotted for the NASA Lewis gear fatigue test apparatus 

A profile modification design method to minimize the vibra- 
tion excitation arising from a pair of meshing gears is re- 
viewed and extended. Families of tooth loading functions re- 
quired for such designs are developed and examined for poten- 
tial excitation of individual tooth vibrations. The profile 
modification design method is applied to a pair of test gears in 
the NASA Lewis fatigue test apparatus using an approximate tooth 
stiffness calculation provided by NASA. 

The above described methods are used to compute the bearing 
support force spectra and temporal histories of the force trans- 
mitted by a pair of teeth and the tooth root stresses for two 
sets of profile modifications — a conventional modification con- 
sisting of linear tip relief, and the above mentioned profile 
modification for minimum vibration excitation. Significant dif- 
ferences in the dynamic forces transmitted by the teeth and In 
the tooth root stresses are predicted for these two cases. 

These differences are discussed and are related to the gear 
system dynamic properties. 

The above computational methods are carried out, for the 
most part, in the frequency domain using the fast Fourier trans- 
form computational algorithm. The required computer programming 
was carried out by Mr. Robert Fabrizio with supplementary pro- 
gramming by Ms. Caroline Buchman and Mr. Edward Campbell. 



THE EXCITATION 


A time varying component of vibratory excitation arises 
from each meshing pinion-gear pair in a gearing system. For 
each gear element — i.e., pinion or gear — the principal source 
of this excitation component is the deviations of the active 
tooth surfaces of that gear element, under loading, from perfect 
involute surfaces with uniform spacing. Thus, this principal 
source of excitation is most naturally described as a displacement 
type of excitation as opposed to a force type of excitation. This 
displacement excitation gives rise to a time varying component in 
the total force transmitted normal to the tooth surfaces which 
is also dependent on the dynamic parameters of the gearing system — 
i.e., on the overall configuration and the mass and interial 
properties of the system. 

In addition, there are time varying tangential forces between 
meshing pairs of teeth that arise from the friction and sliding 
action of the teeth. The component of vibratory excitation caused 
by these frictional forces is generally believed to be small in com- 
parison with the component due to deviations of the loaded teeth 
from perfect involute surfaces. These frictional forces therefore 
will be neglected in this report. 

Gear elements with perfect, rigid, uniformly spaced involute 
teeth transmit exactly uniform angular velocities. Thus, to pre- 
dict the response of a gearing system to deviations of the loaded 
teeth from perfect involute surfaces, we require the deviations of 
the angular positions of individual gear elements under loading 
from the positions of their perfect, rigid, involute counterparts. 

In this report, we shall use the static transmission error [1-3] 
to describe this component of the vibratory excitation, where the 
static transmission error of a meshing gear pair can be defined 
[4,5] as the deviation 60 from linearity of the angular position 0 
of a gear measured as a function of the angular position of the 
gear it meshes with when the gear pair is transmitting a constant 
torque at low enough speed so that inertial effects are negligible. 
Our analysis to follow will use the lineal static transmission error. 


? = R b 60 » (1.1) 


where R b is the base circle radius [6] of the gear whose angular 
transmission error is 60. When the equations of motion of the 
gearing system are written later in this report, we shall see how 
the lineal transmission error C is used to predict the vibratory 
response of points of interest within the system. 

The static transmission error of a pair of meshing gears is 
made up of contributions from each of the two gears. For each 
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gear, there is a contribution from the elastic deformations 
of the teeth and gear body, and a contribution from 
the geometric deviations of the unloaded teeth from perfect uni- 
formly spaced involute surfaces. Furthermore, in References 4 
and 5, we have found it useful to further decompose these geo- 
metric deviations into a mean component and a random component. 

The mean geometric deviation component for a pinion or gear is 
defined as the tooth surface formed by taking the average of all 
tooth surfaces on the pinion or gear under consideration. The 
random component of the geometric deviation of a tooth surface 
is defined as the deviation of that tooth surface from the mean 
tooth surface. Thus, every tooth surface on a pinion or gear 
has the same mean deviation component; however, the random devia- 
tion component generally will differ from one tooth to the next. 

It is shown in References 4 and 5 that the elastic deformations 
and the mean component of the geometric deviations of the teeth 
give rise to the tooth meshing harmonics of vibratory excitation 
and response, whereas the random component of the geometric de- 
viations gives rise to the rotational harmonics, and especially, 
to the so-called sideband components of the spectrum which occur 
at the tooth meshing harmonic frequencies plus and minus one or 
a few rotational harmonic frequencies. The random component of 
the deviations provides no contribution to the tooth meshing 
harmonics. Proofs and full discussion of these facts can be 
found in References 4 and 5- 

In high quality, highly loaded, ground aerospace gearing, the 
random component of the deviations of the tooth faces from perfect, 
uniformly spaced involute surfaces is generally small in compari- 
son with the mean component of the deviations and the component 
due to elastic deformations. Hence, in the analysis to follow, we 
shall mainly concentrate on the mean deviation of the tooth sur- 
faces and the component caused by elastic deformations. A full 
discussion of the random component of the deviations and its 
effects on the excitation spectrum can be found in References 4 
and 5 • 


Fourier Series Representation of 
Tooth-Meshing Harmonic Components of the Excitation 

Computation of the response of vibratory systems to dynamic 
inputs generally is most easily carried out in the frequency 
domain, where temporal convolution or Duhamel integrals are re- 
placed by simple multiplications. Since the excitation components 
of most interest are the mean component of the deviations of the 
tooth surfaces from perfect involute surfaces and the component 
due to elastic deformations, each of which gives rise to a peri- 
odic excitation with period equal to the tooth spacing interval, 
the appropriate description of the static transmission error of 
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these components for later use in the equations of motion is their 
Fourier series representation. To generate the Fourier series 
representation of these static transmission error components, we 
shall use for our independent variable the quantity 


X = R b e 


5 


( 1 . 2 ) 


where 0 is the nominal angular position in radians of one gear of 
the pair whose base circle radius is Rv, . Thus, the same lineal 
variable x describes the nominal positions of both gears in a 
meshing pair, since the product R 0 is the same for both of the 
gears . 

A thorough analysis and discussion of the Fourier series repre- 
sentation of the static transmission error can be found in Ref- 
erences 4 and 5. Thus, in this section, we shall summarize the 
formulas required for computations of the Fourier series represen- 
tation of the mean or deterministic component £ (x) of the static 
transmission error of a generic meshing pair of m spur gears with 
appropriate reference to the derivations of these formulas. We 
shall then provide a brief treatment of the contribution from 
tooth-spacing errors. 

From Eqs. ( 6 0 ) and ( 6 1 ) of Reference 4, we see that when the 
random component C r (x) of the static transmission error is ne- 
glected, we can express the resulting mean component of the static 
transmission error as 


? m (x) ' ? W Cx) . 


+ C (1) (x) 
m / 


+ S (2) (x) 
m 


(1.3) 


where ? w (x) Q is the 
e las ticae format ions 




and 




(loading dependent) component resulting from 
of the teeth and gear bodies, and 

each represents the contribution to the static 


transmission error arising from the mean tooth face deviations on 
gear (•). Numbered superscripts in parentheses designate one gear 
or the other in a meshing pair. Superscript W on the loading 
dependent component represents dependence of that component on force 
W transmitted by the mesh, and subscript 0 designates that the 
transmitted load is assumed to be the constant value W . 


Following Eqs. (64) through (67) of Reference 4, we can ex- 
press the Fourier series representation of the mean component of 
the static transmission error as 


4 


oo 


II 


! m (x) ' ZXn 

n=-a> 


(1.4) 


where, from the linear relationship on the right-hand side of 
Eq. (1.3), the Fourier series coefficients a mn of S m (x) may be 
expressed as 


a 


mn 


A/2 

■ 

S (x) exp(-i2iTnx/A) dx 

-A/2 


= a 


Wn 


+ a 


( 1 ) 

"mn 


+ a 


( 2 ) 


mn 


(1.5) 


and where a , and are, respectively, the Fourier expan- 

sion coefficients of the load-dependent component S T .(x) n and the 

M \ (2) wo 

mean deviation components (x) and '(x) from the teeth of 

gears (1) and (2) of the pair — i.e.. 


a 


Wn 



a 


(*) = 
mn 



C TT (x) exp (-i2Trnx/A) dx 
w o 

£^*^(x) exp(-i2imx/A) dx . 
m 


( 1 . 6 ) 


(1.7) 


Let us concentrate now on determining a useful set of expres- 
sions for computing the Fourier series coefficients a -^ n of the 
loading dependent component of the static transmission error . 

From Eq. (68) of Reference 4, we see that can be expressed as 


“Wn = W ,“(l/K)n> n = 0, ± 1, ± 2, ••• (1.8) 

where the constant force W Q transmitted by the mesh is defined in 
the direction determined by the intersection of the base plane and 
axial plane, and are the Fourier series coefficients of the 

reciprocal of the total mesh stiffness. According to Eqs. (77) 
and (. 87 ) of Reference these latter Fourier series coefficients 
can be computed to any desired degree of accuracy from the expres- 
sion 


5 



1 M /M+l\ 
( 1/K)n £?o l £+1' 



[£] 


3 


(1.9) 


where the accuracy of this expression Increases with Increasing M. 
According to Table I on p. 1761 of Reference 5, M should be taken 

at least equal to M = 3 for accurate calculations of a,, for 

/M+l\ (1/K)n 

spur gears. The quantities f are the binomial coefficients 


£+1 



(M+l ) ! 

(£+1)1 (M-£ ) ! 9 


( 1 . 10 ) 


whereas, in Eq. (1.9), we have also used the definition 



where & n is Kronecker's delta 
n,u 

i = j 


otherwise, (1.12) 

and where the vertically centered asterisks in Eq. (1.11) denote 
discrete convolutions defined by the right-hand equality in 
Eq. (G12) on p. 1428 of Reference 4. Notice that for £ = 1,2,3, .. 
we have 
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f 


The quantities are the Fourier series coefficients of the 

total mesh stiffness a§ defined by Eq. (75) of Reference 4 for 

SL-1, and K T is the mean value of the total mesh stiffness as defined 

by Eq. (F.2) on p. 1427 of Reference 4. These latter two quantities 
can be computed from the spur gear local tooth pair stiffness K s (z) 

which is described in detail in Appendix A of this report. Since we 
are considering here only spur gears, which have zero helix angles, 
we have in these cases from Eq. ( C 4 ) of Reference 4, (L/A) = 0. 

For spur gears, we therefore have from Eq. (86) of Reference 4, 


a Kn _ K TC [0,(nL/DA)] 

k t k tc (o,o) 


(1.14) 


whereas, from Eq. (85) of Ref. 4, we have 


k t = K tc (0,0) , (1.15) 

which reduces our problem to determination of K TC [0 , (nL/DA )] . Let 
us define the spur gear local tooth pair stiffness K s (z) as 
F/2 

Kg(z) 4 | K TC (y,z)dy , (1.16) 

-F/2 

where F is the face width of the gear and K T ^(y,z) is the local 

tooth pair stiffness per unit length of line of contact expressed 
in the tooth coordinates defined by Eqs . (15), (16), and (19) of 
Reference 4. From Eq. (82) of Reference 4 and Eq. (1.16) above, 
we then have 


K TC [0, (nL/DA)] 


and 


D/2 

I 

-D/2 


K s (z)e- i2,,(nL/DA)z dz 


D/2 

K tc (0,0) = j K s (z)dz, 
-D/2 


(1.17) 


(1.18) 


where D is the total height of the active tooth surface as illus- 
trated in Fig. 4 of Reference 4. Coordinate z is defined in 
Appendix A, and L is the length of line of contact illustrated in 
Fig. 2 of Reference 4 which is related to D by 
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L 

D 


CSC(f> , 


(1.19) 


according to Eq. (C8) on p. 1426 of Reference 4 where <p is the 
pressure angle. According to Eqs . (A5) and (A6) on p. 1424 of 
Reference 4, A is the base pitch which is related to the circular 
pitch A c by 

A = A c (R b /R) = A c cos<|> (1.20) 

where R, is the base circle radius and R is the pitch circle 
radius. The overall parameter L/DA that appears repeatedly above 
therefore can be expressed in terms of the pressure angle <j> and the 
circular pitch A q by 


E csciji ctn<f> ( -\ o I 

DA A * ^ ; 

c 

Once a value of M CM^3) -is chosen, Eqs. (1.8) through (1.21) deter- 
mine the Fourier series ooeffioients a ^ n > n=0, ±1, ±2, • * * of the 

load dependent component of the static transmission error in terms 
of the loading W Q carried by the mesh , the local tooth pair stiffness 

Kg(z) over the range -D/2 z ^ D/2, and the design parameters <J> 

and A . 
c 


We turn now to determining an analogous set of expressions 

( • ) 

useful for computing the Fourier series coefficients a of the 

mean deviation components of the static transmission error for 
either gear (•) = (1) or (2) of a meshing pair. According to 
Eq. (94) of Reference 4, the mean deviation Fourier series coeffi- 
cients can be expressed as 


(•) ... 


.(•) 


a' ' =a,, /T .>. *a ' ' , n= 0. ±1, ±2, 

mn ( 1/K)n m , n s s * * 


( 1 . 22 ) 


where the coefficients are given by Eqs. (1.9) through 

(1.21), and where the vertically centered asterisk again denotes 

the discrete convolution defined by the right-hand equality in 

Eq. (G12) on p. 1248 of Reference 4. The coefficients 

( • ) . . 
a are determined from the mean profile modification of gear (•) 

m'n 

by first using Eq. (92) of Reference 4 applied to the case of spur 
gears [where = 0 and (L/A) = 0] — i.e. , 

a ( ’ } = ^ [0, (nL/DA)]. (1.23) 

m’n L ' 
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Since our interest here is spur gearing, we shall now restrict the 
general theory presented in References A and 5 to the case where the 

mean tooth face modification n/ ^(y,z) of gear (•) (defined on 
p. 1417 of Reference 4) is independent of the axial coordinate y 
illustrated in Pig. 4 of Reference 4. That is, let us define for 
spur gears the mean profile modification of gear ( • ) as 


m 


O ’ ^ ( Z ) = it/ ^ ( • , Z ) 


(1.24) 


where ir/ ^(*,z) is the tooth face modification rr/ ^(y,z) that we here 

assume to be independent of axial location y. Let us further define 
for spur gears the stiffness weighted mean profile modification 
(•) 


mi. '(z) as the axial integral of the stiffness weighted tooth face 

IA.O / m \ 

modification m^.^ (y,z) expressed by Eq. (48) of Reference 4 — i.e.. 


m KS >(z) = 


P/2 

’ ( * ) 

m Kc ’ z ^ d y 

-F/2 
P/2 


K TC (y,z)m^’ ^ ( • , z ) dy 


-P/2 


= Kg ( z )nig ^ (z) , 


(1.25) 


where the last line follows from Eqs . (1.16) and (1.24). Setting 
g = 0 in Eq. (91) of Reference 4, we see that 

Y N 

in^c ^ [0 , (nL/DA ) ] can b'e expressed using Eq. (1.25) as 


^KC } E 0 j ( nL/DA ) ] 


' D/2 K s ( z )m< • 5 ( z Je- 1211 (nL/D& > z dz . 
-D/2 


(1.26) 


Equations (1.22) y (1.2Z) 3 and (1.26) determine the Fourier series 

. . ( • ) 

coefficients af , n = 0, ±1, ±2, ••• of the contribution to the 
mn 3 3 J 

static transmission error from the mean -profile modification of 
gear (•), (•) = (1) or (2), in terms of the coefficients 
n = 0, ±1, ±2, *•* determined by Eq. (1.9) 3 the local tooth pair 
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stiffness K q (z) described by Eq. (1.16) 3 the mean profile modification 
nig (z) of gear (•), and the design parameters given by Eq. (1.21). 

When the Fourier series coefficients of the load-dependent 
component of the static transmission error , Eq. (1.8 ) , are added 
to the Fourier series coefficients of the mean profile modification 
component from each gear 3 Eq. (1.22) 3 (•) = (1) or (2) 3 as shown 

in Eq. (1.5) 3 we have a complete expression for the Fourier series 
coefficients of the tooth meshing harmonic components of the static 
transmission error of the meshing gear pair. These formulas, and 
the ones to follow, take into full account the alternating numbers 
of teeth in contact as the gears rotate. 

Fourier Series Coefficients of Rotational Harmonic 
Components Caused by Tooth-Spacing Errors 


A thorough discussion of the contributions of manufacturing 
errors that cause variations in the running surfaces of gear teeth 
from one tooth to the next can be found in References 4 and 5- 
That treatment is capable of predicting the rotational harmonic 
contributions caused by any errors of this type. However, in high 
quality ground aerospace spur gears, we would expect the dominant 
manufacturing error component to be that associated with tooth spac- 
ing errors. Formulas for the rotational harmonic contributions 
caused by tooth spacing errors are summarized below. 


The tooth spacing errors on gear (•) of a meshing pair, 

(•) = (1) or (2), give rise to a periodic contribution to the static 
transmission error, where the period is the base circle circum- 
ference of the gear. If N' ' denotes the number of teeth on the 

gear and A the base pitch, then that period is N' A. Thus, there 
are two sets of such rotational harmonics generated from each pair 
of meshing gears, one set from each gear, unless the numbers of 
teeth on the two gears is the same. 


Let e 


(•) 


denote the accumulated tooth spacing error of tooth 


j of gear (•) as defined by Eq. (107) of Reference 5 with 

(•) 


E <-> = b (-) 

J - b j,oo- 


Let B (n) denote the finite discrete Fourier trans- 


form of the sequence e : 

J 


(•) 


, j = 0,1,* • * ,N-1 : 


B (#) (n) = 


N ^ * ^-l 


N 


1 


i 

j=0 


(.) -i2irnJ/N (,) , n=0 , ±1 , ±2 , • • • . (1.27) 


Then, we see from Eqs . (52) and (53) of Reference 5 that the Fourier 

( • ) 

series coefficients of the component of the static transmission 
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error caused by the tooth spacing errors of gear (•) can be 
expressed by the simple multiplicative relationship 


a'; ) = | 

rn o o 




, n = 0,±1,±2, 


( 1 . 28 ) 


where, here, n denotes the rotational harmonic of gear (•) where 

7 .) 

the period of the fundamental component n = 1 is ' A. Prom 
Eq. (110) of Reference 5 , we see that the mesh transfer function 

for tooth spacing errors $ ( vvy ) appearing in Eq. (1.28) can 

be expressed as 00 'lA 


0 °\? * >A 


n \ _ sin[ (n/N^ * ^ )ttL/A] 


(n/N ( ’ 5 )ttL/A 


V 


w 


0 0 


/ n A sin (n 1 ttL/A ) s in|[(n/N^ A-n’]TrL/A[ 

\N " ^ ’ / n’TTL/A [ (n/N^ * ^ )-n' ]ttL/A 


except 
n 1 =0 


(1.29) 

where we have used g = n/lA ^ and the fact that j fl (x) = sinx/x 
as indicated by Eq. (126) of Reference 5- The function w in 
Eq . (1.29) is defined by Eqs . (93) and (9*0 of Reference 5 — i.e.. 


w 


o o 


■ -•) 


= 2 rect 


1 / n 


2 ^ 


- n, )]HK#r-.n')]j 


(1.30) 


where from Eq . (90) of Reference 5, we see that 

1 , 


rect[x] = 


x i < \ 


(1.31) 


0 , 


x | 5- ~ 


In the above formulas, the quantity L/A is the transverse contact 
ratio which, from Eq. (148) of Reference 5, we see can be expressed 
in terms of fundamental design parameters by 


L 

A 


DN ^ ^ csccj) 


2ttR 




(1.32) 


ii 


H 

(•) (•) 

where, as before, D, N , <j>, and R£ are, respectively, the 

active depth of the teeth, the number of teeth on gear (•), the 
pressure angle, and the base circle radius of gear (•). 

Use of the above formulas is illustrated on pp. 1774 through 
1780 of Reference 5* A discussion of the general properties of 
the excitation spectra of spur and helical gears may be found on 
pp. 1781 through 1785 of that same reference. 
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THE EQUATIONS OF MOTION 


We turn now to writing the system equations of motion in a 
form useful for computing the desired response metrics from the 
excitation characterization described in the preceding section. 

The system of Interest is the "four-square" gear fatigue test 
apparatus located at the NASA Lewis Research Center. A schematic 
diagram of this system is shown in Pig. 1. The system contains two 
test gears and two slave gears as shown. The lightweight test gears 
are connected to the relatively massive slave gears by short very 
stiff shafts. Torque is applied between the shafting and one of 
the slave gears using a pressurized-oil torque-applier . There is 
no accumulator associated with the torque applier. The system is 
firmly supported. 

A displacement type of excitation takes place at each of the 
two meshes I and II. Because of the relatively large stiffnesses 
of shafting, gear bodies, and bearing support structure we shall 
assume that the (displacement) excitation at the meshing points 
is "taken up" by tooth elastic deformations and relative motion 
between the shafting and slave gear at the fluid loaded torque 
applier. That is, shafting, gear bodies , and bearing supports shall 
be modeled as rigid members with gear bodies and shafting possessing 
inertia. Gear teeth shall be modeled as elastic members. The 
torque applier shall be assumed to apply constant torque x plus a 

dissipation torque c (0-0,) proportional to the angular velocity 

difference between the slave gear and its shafting. Thus, the 
system is assumed to have three degrees of freedom, 0 , 0 2 , and 0 3 

as shown in Fig. 1. Since some vibratory energy associated with 
motions 0 3 and 0 2 is necessarily dissipated in the bearings, 

viscous damping terms are directly associated with these two degrees 
of freedom. For completeness, a viscous damping term also will 
be associated with 9 . 

3 

Because bearing supports are assumed to be rigid, each gear 
must move in pure (generally unsteady) rotation. Thus, the instan- 
taneous forces Wj and W^j transmitted by meshes I and II must have 

equal and opposite reactions at the bearing supports as illustrated 
by the test gear supports in Fig. 1. Since we are dealing with 
gears of nominal involute design, the directions of mesh forces 

and and their bearing reactions, remain fixed and parallel to 

their respective planes of contact. However, the magnitudes of 
these forces generally will possess temporal variations. 

All gears in the system have the same base circle radius R, 
which is the radius Instrumental in generation of the gear torques 
by forces Wj and W^-j. I designates the moment of inertia of the 

left-hand slave/test gear pair and shafting whose angular position 
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Torque Applied 
to Geer by 
Shaft * t 

Torque Applied 
to Shaft by 

r-T.+ c.W,-*,) 


FIG. I. SCHEMATIC OF NASA LEWIS GEAR FATIGUE TEST 
APPARATUS. 
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I 


is e . I designates the moment of inertia of the test gear and 
1 2 

shafting whose angular position is 0 , and 1^ designates the moment 
of inertia of the slave gear whose angular position is 6 . Dissi- 
pation (dashpot) constants c , c^, and c ^ are associated with 
0,0, and 6 respectively. 

12 3 

Thus, the differential equations of motion of the three 
rotational elements are 


I 0 

l l 

+ 

c 0 i 
i 

= RW i: j- 

- RWj = R(W II -W I ) 

(2.1) 

I 0 

+ 

c 6 

= -rw t 

+ x + c ( 0 -0 ) 

(2.2) 

2 2 


2 2 

I 

a a 3 2 


I 0 

+ 

c $ 

= rw tt 

- x - c ( 0 -0 ) . 

(2.3) 

3 3 


3 3 

II 

a a 3 2 



The moments of inertia 1,1, and I , radii R, applied torque x , 

123 a 

and dissipation constants c lS c 2 , c 3 , and c all are assumed to 

be independent of time, whereas, angular positions 0 , 6^, and 0^ 

and forces and are assumed to be generally time dependent. 

Angular positions 0 , 0^, and 0 are the vibratory components of 

the true angular positions of the gears. 

Let us now relate the instantaneous mesh forces and to 

the transmission errors of the appropriate gear pairs and the 
angular positions 0 of the gears. The transmission error is mea- 
sured in the axial plane in a direction parallel to the plane of 
contact of the gear teeth, and is defined as positive when it is 
"equivalent" to removal of material from the tooth surfaces 
[p. 1410 of Reference 4]. Thus, if we denote the (dynamic) trans- 
mission error of meshes I and II by ^ and respectively, then 

for the teeth to be in continuous contact we require the following 
two compatibility relations to be satisfied [p. 1410 of Reference 4] 
as may be seen from Pig. 1: 

R0 + R0 = 4 t (2.4) 

1 2 1 

-R0 - R0 = C TT . (2.5) 

Let ? denote the (dynamic) transmission error or 

associated with either of the two meshes illustrated in Fig. 1, and 
let C|j, denote the total compliance of the teeth in the same mesh. 

C|p varies periodically with time due to variations in tooth-pair 

stiffness and the numbers of teeth in contact. Let W denote the 
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force transmitted by the mesh which also can vary with time. 

Then can be expressed [p. 1410 of Reference 4] as 

= C T W (2.6) 

where and Z each denotes the component of the transmis- 

sion error from one gear in the mesh due to geometric devia- 
tions of the active tooth surfaces from perfect involute sur- 
faces [p. 1410 of Reference 4], The average value of and 

(2) 

Z , where each value is an average over all teeth on the appro- 
priate gear (1) or (2), yields the mean transmission error com- 
ponents and contained earlier in Eq. (1.3) and de- 

scribed in Reference 4. 

To enable us to use a transfer function approach to effect 
a solution to the gear system dynamic problem, let us now de- 
compose C T and W into their time average components C T and W , 
and their fluctuating components <50^ and 5W — i.e.. 


C T = c T + 

sc T 

(2.7) 

w = w + 
0 

6W . 

(2.8) 

We then have 



c T w = c T w 

+ 6C T W + <5C t 6W 

(2.9a) 

* V 

+ w sc 
o 1 

(2.9b) 


where, in the second line, we have neglected the product 

of fluctuating quantities. If we now define a loading dependent 
component of the static transmission error as 


Cw 0 ^ W ° 6C T » (2.10) 

we see that since (1/K) = C T , this definition is the same as that 
of £ w (x) q in Eq. (1.3) whose Fourier coefficients are given by 
Eq. (1.8) — except that the dc component W 0 C,p of ?^(x) Q de- 
scribed by the term n = 0 in Eq. (1.8) is not included in the 
definition, Eq. (2.10). Using the approximation given by 
Eq. (2.9)-, we may now combine Eqs. (2.6), (2.9), and (2.10) to 
give 
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C D = C T W + c; + c (1) + e (2) 
0 


= C T W + c’ , 


( 2 . 11 ) 


where In the second line we have defined 







( 2 . 12 ) 


which is the static transmission error discussed in the first 
section of this report except for the dc component W o C T not 

included in Eq. (2.12). Equation (2.11) provides an approximate 
relationship between the dynamic and static transmission errors 
that would be exact if the product <SC T SW were added to its 
right-hand side. 

Using the approximation given by Eq . (2.11), we therefore 
can express the dynamic transmission errors of each of the two 
meshes illustrated in Fig. 1 by 


and 


- c lWj + c’ 


(2.13) 

(2.14) 


where and denote respectively the time-average values of 

the compliances Cj and Cjj of meshes I and II, Wj and denote 

the time-varying forces transmitted by these meshes, and Cj and 
Cjj denote the static transmission errors of these meshes as 


defined by Eq. 
Eqs . (2.5) and 

(2.12). 
(2.14) , 

Combining Eqs. (2.4) and (2.13), 
gives 

and then 

R(0 +e ) = 
1 2 

c T Wi + 

4 

(2.15) 

and 




-R(0 +0 ) = 

1 3 

°II W II 

+ *ii > 

( 2 . 16 ) 

which we can solve algebraically for Wj and to give 


W T = (C T )“ 

1 [R(e +e 

)-?;] 

(2.17) 
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and 


W n = “ Cc ii ) “ 1 C R( - e 1 +e 3 )+? ii ] * (2.18) 

Mesh forces W T and W TT , system displacements 0,0, and 0 , and 

12 3 

transmission errors and Cj-j- in Eqs . (2.17) and (2.18) are all 

assumed to be time-varying quantities. Finally, by substituting 
the expressions for W^. and given by Eqs. (2.17) and (2.18) 

into Eqs. (2.1) through (2.3), we obtain after minor rearrangements 


T 9 +c 
i i 


0 +R 2 [(C t )" 1 +(C tt ) _! ]9 +R 2 (C T ) ” 1 0 
1 1 l-J- 1 -1-2 


+R 2 (C ii ) _ 1 0 3 = R(C I )" 1 <; ; j.-R(C II )“ 1 ^ I 


(2.19) 


r 2 (c t ) _1 0 +i e +(c +c )e +r 2 ( c t ) - 1 e -c e 

I 122 2a2 1 2a3 


■ V R(c i r ‘ ? i 


R 2 (c TT )“‘e -c e +i e +(c +c ) 0 +r 2 (o tt ) _1 0 

ii i a 2 33 3 a 3 II 3 


T a R(c n ) 1? n • 


( 2 . 20 ) 


( 2 . 21 ) 


Equations (2.19) through (2.21) are the desired differential 
equations of motion for the system (displacement) responses 
0,0, and 0 with the (time-varying) static transmission errors 

12 3 j 

SjL and of meshes I and II as the excitation. C-j- and are 

defined by Eqs. (2.12) and (2.10), and may be written explicitly as 


(,T a /R)6C I + 




(2.22) 


and 



(V R)6C II 


+ C 


( 1 ) 

II 


.(3) 

■II 


(2.23) 


where t /R is the time-average force W Q transmitted by the two meshes, 

6C t and <5C TT are, respectively, the fluctuating components 6C m of 
l 11 (1) (2) 1 

the compliances of meshes I and II, and are the geometric 
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components of the transmission errors of gears (1) and (2) of 
mesh I, and i^j^and are the geometric components of the 


transmission errors of gears (1) and (3) of mesh II. See Pig. 1. 
The derivation of Eqs . (2.19) through (2.21) illustrates a speci- 

fic application of the general methodology described on p. 1410 
of Reference 4. 


Solution Using Transfer Functions 


We shall now solve the set of Equations (2.17) through (2.21) 
to yield algebraic expressions for the Fourier transforms of the 
time-dependent mesh forces W^- and W-^-j- in terms of the Fourier 

transforms of the static transmission errors and (or, 

equivalently, their Fourier series coefficients). We shall begin 
by taking the Fourier transforms of each of the five Equations 
(2.17) through (2.21). In general, we denote the Fourier trans- 
form of a variable by placing caret over that variable — e.g., for 
0 (t), 0 (t), or 0 (t), we define 

12 3 


0 (to ) ^ 


0 (t )e 


-irnt 


dt 


(2.24) 


for W^(t) and w jj(t), we define 


W(m) ^ | W(t)e dt , 

— oo 

for Cj(t) and Sjj(t), we define 

OO 

8' (co) & | ? '(t)e- lwt dt , 


(2.25) 


( 2 . 26 ) 


and, finally, for our applied dc torque t , we define 

cl 


V u) = { 


x e lwt dt. 
a 


(2.27) 


The Fourier mates to Eqs. (2.24) through (2.27) are 


8(t) = ! 


6(m)e ia)t dm 


J 

—CO 


( 2 . 28 ) 


W(t) = 


2tt 


W(w)e^ wt dm 


(2.29) 
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and 


C'(t) = 


' a 2 it 


2tt J 


$ i , > icot , 

£ (co) e dco 


~ > \ icot , 
x (co) e dco , 

B. 


( 2 . 30 ) 


(2.31) 


where, from this last relation, we may conclude that 

T a (oi) = 2 7 T t a 6 ( u ) (2.32) 

since x is a constant torque. In Eq . (2.32), <$(co) denotes the 

cl 

Dirac delta function, which is zero everywhere except at to = 0, 
and whose integral over co is unity. Prom Eq. (2.28), we may also 
conclude that 

§<t) • h [ 

and 

« (t) - h [ 

from which it follows that ico9(co) and -co 2 0(co) are, respectively, 
the Fourier transforms of 0(t) and 9(t). 


ico0 (coje^ 00 ^ dco 


-co 2 0(co)e 1 dco. 


(2.33) 


(2.34) 


Utilizing Eqs . (2.24), (2.26), (2.27), (2.33), and (2.34), we 
may now write the Fourier transforms of Eqs. (2.19) through (2.21) 
as 


i-I co 2 +ic co+R 2 [(C t )" 1 + (C tt )“ 1 ] l© +R 2 (C t ) -1 0 
(.1 1 1 11 J 1 1 2 

+ r 2 (c ii )~ 1 0 3 = r(c i )“ 1 8j-r(c ii ) -1 ?j I 

R 2 (C t ) -1 0 +[-I co 2 +i(c+c )co+R 2 (C T ) _1 ]§ 

1 1 2 2 3. 1 2 

- ic a“ S , - + 

and 

R 2 (C tt )“ 1 9 - ic co 9 +[-I co 2 +i ( c +c )co+R 2 ( C TT ) ” 1 ] 9 

ix 1 a 2 3 30- xx 3 


(2.35) 


(2.36) 


(2.37) 
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w 


Equations (2.35) through (2.37) are a set of linear algebraic 

equations for the Fourier transforms of the responses 6,6, ^tnd 

6 in terms of the Fourier transforms of the excitation^ tJ t, 
s /n a -L 


3 

and 


c' . Let us define H. .(co) as the response 0. due to unit 

11 At 1 
excitation = 1, where i = 1,2, or 3 and j = I or II. Then 

for io ^ 0, the set of Equations (2.35) through (2.37) can be 
written in matrix form. 


[A ±j . ( w)] [H ±J (co) ] = [U ± j ( w ) ] , “ * 0 


(2.38) 


where the condition id ^ 0 arises from the fact that x (oj) has a 

cl 

component only at io = 0. From Eqs. (2.35) through (2.37)» we 
see that [A^tco)] is a symmetric square matrix with elements 


A 

l l 

= 

3 

H 

H 

1 

A 

1 2 

= 

R 2 (C. 

A 

1 3 

= 

R 2 (C 

A 

2 1 

= 

A 

1 2 

A2 2 

= 

-I uo 
2 

A 

2 3 

= 

-lc a l 

A 

3 1 

= 

A 

1 3 

A 

3 2 

= 

A 

2 3 

A 

3 3 

= 

-I U) 
3 


'II' 


, - i 


, — i 


!!■ 


L“ 1 


L(c +c )w+R 2 (C tt ) *, 

3 a- -1 -L 

whereas the excitation matrix [U. . (u))] can be written as 

1 3 


(2.39) 


C TJ 1 j( f °)] = 


R(C I )" 1 

R(C I ) _1 


-R(Cn )- 1 


-R(C II )" 1 


(2.40) 
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We wish now to express the Fourier transforms of the mesh 
forces W-j- and W^ in terms of the solutions of the matrix 

Eq. (2.38). Forming the Fourier transforms of Eqs . (2.17) and 
(2.18) and using the definitions, Eqs. (2.24) through (2.26), we 
have 


and 


Wj = (c I )' 


W 


H 


. r a a a f n 

[^vv-'i] 


- - (o n ) ' i [ R( t + s )+ ?d- 


(2.41) 


(2.42) 


Recognizing that H 1 - J -(m) and H^Coo) are the responses 6^(00) and 
to a unit input C^(w) = 1 at mesh I, and that and 


e (u) 

2 

H ott (w) are the responses 0 (00) and 0 (w) to a unit input 

a ^ -L X 1 2 

Cjj(w) = 1 at mesh II, it follows from the linearity of the system 

of Eqs. (2.35) through (2.37) and the definitions of H^-(m), 
i = 1 and 2 and j = I and II that from Eq. (2.4l) we have for the 
Fourier transform of W^(t) , 




= (c I )" 1 


^{R[H 1I (.o3)+H 2I (oj)]-l}6j 


+ ( w ) + H 2 (w) ] 

Similarly, we may express W^-j- as 


a) / 0 


(2.43) 


W. 


II 


= -(C II )" 1 |R[H 1I (a))+H 3I (m)]8-[ 

+ {R[H 1 II (.m)+H 3 II (a))]+l}Cj I ) , 


w ^ 0, 


(2.44) 


Equations (2.43) and (2.44) are valid for all w except to = 0. 
Since our central interest is the force history in mesh I, we 

A 

shall evaluate from here on only the Fourier transform of that 

force history. From Eq. (2.43), we see that once H^-j-(to), H 2 j(to), 
Hm(co), and H 2 jj(to) are evaluated, the computation of Wj(co) from 
C-j-(w) and Cjj(co) is a very simple matter. 


Evaluation of Transfer Functions 


Let us now proceed to evaluate the transfer functions H^.tio), 

i = 1 and 2, j = I and II, which are solutions to the matrix equa- 
tions defined by Eqs. (2.38) through (2.40). Using Cramer's rule 
[e.g. , p. 12 of Reference 7], we recognize that any of the required 
Hij (w) can be expressed as 


22 


> 


(2.45) 


N (03) 

■ Tfer 


where D(oo) is the determinant of the coefficient matrix [A^tco)] 

whose elements are given by Eq. (2.39). The determinant D(ai) is 
readily evaluated to give 


D(o>) = \A ± . 


A 


A 

A 

l 

l 

1 2 

1 3 

A 


A 

A 

2 

l 

2 2 

2 3 

A 


A 

A 

3 

l 

3 2 

3 3 

A 

(A 

A -A 

A ) 

l l 

22 33 23 

3 2 

-A 

(A 

< 

i 

< 

A 

2 

i 

1 2 3 3 1 

3 3 2 

+A 

(A 

A -A 

A 

3 

i 

12 23 13 22 

A 

(A 

> 

1 

> 

to 

) 

i l 

22 33 23 


-A 

(A 

A -A 

A 

i 

2 

1 2 3 3 1 

3 2 3 


(2.46a) 


(2.46b) 


+A (A A -A A ) , 

13 1223 1322 


(2.46c) 


where we have used the symmetry property A.. = A . in writing this 
last expression. Using Cramer’s rule, it follows from Eqs . (2.38) 


N 1I (w) = 


R(C I )“ 1 

R(C I )" 1 


11 

(to) as 

1 2 

A 

1 3 

2 2 

A 

2 3 

3 2 

A 

3 3 


(2.47a) 


= R(C t ) _1 [(A A -A A )-(A A -A A )] 

1 2233 2332 1233 1332 


= R(C t ) -1 [A (A -A )+A (A -A )] (2.47b) 

I 33 22 12 32 13 23 

= R(C t ) _1 [A (A -A )+A (A -A )], (2.47c) 


33 22 12 


2 3 1 3 2 3 


where we have again used A 
larly, we have 


ij 


= A . ^ in writing the last line. Simi- 
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N 2I (co) A 


(2.48a) 


1 1 

R(C I ) _1 

A 

1 3 

2 1 

R(C I )- 1 

A 

2 3 

3 1 

0 

A 

3 3 


= R(C t )" 1 [A (A -A )+A (A -A )] 

r 33 11 21 31 23 13 


T 

'I' 


= R(C t )" 1 [A (a -a )+A (a -A )], 

1 33 11 12 13 23 13 


N 1II (^ ) 


-r(c tt )- 1 a 
11 1 2 


1 3 


A 


2 2 


2 3 


-R(C T )" 1 A 
II 3 2 


A 


3 3 


and 


r(c tt )- 1 [a (A -A )+A (A -A )] 

22 13 33 23 32 12 


II' 
•II : 


R(C tt )- 1 [A (a -a )+A (a -a )], 

11 22 13 33 23 23 12 


N 2 h (w) 


A -R(C tt ) ” 1 

ii ii 

A 0 

2 1 


Ai 3 

A 


2 3 


A 


3 3 


R(C TT )- 1 [A (A -A )+A (A -A )] 

11 21 33 13 23 11 31 J 

R(C tt )- 1 [A (A -A )+A (A -A )] 

11 12 33 13 23 11 13 


(2.48b) 

(2.48c) 

(2.49a) 


(2.49b) 

(2.49c) 


(2.50a) 

(2.50b) 

(2.50c) 


The symmetry property A^. = A ± has not been used in the "a and b" 

versions of Eqs. (2.46) through (2.50); however, this symmetry 
property has been used in the "c" version of Eqs. (2.46) through 
(2.50). 


From Eq. (2.43), we see that we require the sums 
Hu ( w ) + H 2i (c °) and H^jj(io) + ^ 211 ^°^* According to Eq.(2.45) 
can express these sums as 
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3 


we 



I 


i* 


H i;e (w) + H 2 I (co) = 


N 1I (w)+N 2I (w) 

dTO 


and 


H 1 II (w) + H 2 II (to) = 


N i;ri (u))+N 2II (a») 


(2.51) 


(2.52) 


Prom Eqs. (2.47) and (2.48), It follows that the numerator In 
Eq. (2.51) can be expressed as 


N 1I + N 2I = R(0i)' , [A,, ( A, 1 +A 2i -A 1I -A !1 ) 

+ (A -A HA -A )] (2.53a) 

1 3 2 3 3 2 3 1 

= R(C t )- 1 [A (A +A -2A ) -( A -A ) 2 ], (2.53b) 

1 33 11 22 12 13 23 


where A. . = A., has been used In the second line, but not In the 
a J <] a 

first line. In an analogous manner, we have from Eqs. (2.49) and 

(2.50) 


n 1 tt +n 9tt = R(C tt ) -1 [A (A +A -A -A ) 
111 211 II 23 11 32 12 31 

- (A -A ) (A -A )] 

3 3 1 3 2 2 2 1 


(2.54a) 


= R(C tt ) _1 [A (A +A -A -A ) 

11 23 11 23 12 13 

- (A -A ) (A -A )], 

3 3 1 3 2 2 1 2 


(2.54b) 


where A^. - A^ ^ has again been used in obtaining Eq. (2.54b) from 
Eq. (2.54a). 


Equations (2.46c), (2.53b), and (2.54b) are the relations 
required for evaluation of the sums H^j(ui) + H 2 j(w) an< ^ 

Hin(co) + H 2 jj(w), using Eqs. (2.51) and (2.52), for substitution 

into Eq. (2.43). The coefficients A., are given by Eq. (2.39). 

ajJ 

The results of these straightforward but tedious algebraic mani- 
pulations are 

N 1 I (oi) + N 2 I (m) = R(C I )- 1 |(I i +I 2 )I 3 to 4 

- [(I +1 +1 )R 2 (C TT ) _1 + (c +c +c ) c +(c +c )cja ) 2 

1 2 3 11 1 2 3 a 12^ 

+ i{-[(I +1 )(c +c ) +1 (c +c +c )]w 3 

12 30 . 3120. 

+ (c +c +c )R 2 (C TT )~ 1 ai}) , (2.55) 

123 11 / 
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N^toO + N 2 II (oj) = R(C ii )- 1 {-I 2 I 3 ^ 

+ [(c +c +c )c +c c ]w 2 
1 2 3 a. 2 3 

+ i[ ( I +1 +1 ) c +1 e +1 c ]o) 3 l , (2.56) 

1 2 3 3.23 32 J 

and 

D(w) = -I I I w 6 +{l I R 2 (C tt )' 1 +I I R 2 (C t )" 1 
123 12 II 13 I 

+1 I R 2 [(C t )" 1 +(C tt ) _1 ]+I (C C +C C +C C ) 

23 I II i 23 2 a 3 a 

+1 c (c +c )+I c (c +c ) }o) 4 
213a 312a 

-{(I +1 +1 )R 4 (C T ) _1 (C TT ) -1 +c (c +c ) R 2 ( C TT ) “ 1 

123 I II 1 2 a ±± 

+c (c +c )R 2 (C t ) _1 +[(c c +c c +c c ) ]R 2 [ (C t ) " 1 + (C tt ) “ 1 ] }oo 2 
13a i 23 2 a 3 a l n 

+i ( [I I (c +c )+I I (c +c )+I I c ]o) 5 

'12 3a 132a 231 

-{I C(c +c Q )R 2 (C TT )“ 1 +(c +c a )R 2 (C T )' 1 ] 

+I 2 [c i R 2 (C II )" 1 + (c 3 +c a )R 2 (Cj 1 +c‘J-)3 
+1 [c R 2 (C t )- 1 + (c +c_)R 2 (C t “ 1 +C ~l )] 

3 1 I 2 3 1 11 

+c (c c +c c +c c,)}io 3 

1 2 3 2 a 3 a 

+[(c +c +c )R 4 CC t )" 1 (C tt )- 1 ]w) . ( 2 . 57 ) 

123 I II / 

Low and High Frequency Asymptotic Behavior 


Returning to the expression for W_ given by Eq . 
us define 

(2.43), let 

Qj j(w) A R[H i:[ (u)+H 2I (u)]-1 

(2.58) 

and 


‘ H[H in (»)M 2II (,)] 

(2.59) 

which enable us to express as 


W T = (C i )“ 1 [Q ii (oj)^+Q i u (co) C jp] , 00 * 0 . 

(2.60) 
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Since CyWj, Zj> and is each the Fourier transform of a dis- 
placement, we see from Eq. (.2.60) that Q^Coi) and Q j is 

each a dimensionless transfer function. In general, Q^j (w) 
represents the contribution to the complex displacement amplitude 

C, W. from a unit sinusoidal transmission error amplitude Z '. = 1. 

-*■ i J 

An equation comparable to Eq. (2.60) can, of course, be written to 
replace and interpret Eq. (2.44). 


It Is instructive to examine the low and high frequency 
asymptotic behavior of the two transfer functions Q T (w) and 
Ql jj( u ) in e Q* (2.60). Let us consider first the low frequency 

behavior of Q^^co). Substituting Eq. (2.51) into Eq. (2.58) gives 

after minor rearrangement 

R[N 1 T (w)+N„ T (w)]-D(w) 

- — B® ' < 2 ' 6l > 

Retaining terms proportional to ui and oj 2 in N-^(u>) + ^^(oj) and in 

D(oj) In the numerator of Eq. (2.61), and retaining the term propor- 
tional to oj in the denominator D(oj) of Eq. (2.61), we find that 
the terms proportional to oj in the numerator vanish identically 
leaving for our final asymptotic result as oj -*■ 0, 


Q II (w) 


c c +c c 

OJ 

, 1 2 2 3 

a c +c +c 

12 3 

r 2 (c i )- 1 


oj 0 


( 2 . 62 ) 


which gives the low frequency asymptotic behavior of the transfer 
function Qj-j-(oj) . Notice from Eq. (2.62) that Qjj(^) -> 0 as 

oj 0 , and that the strength of every term in Eq. (2.62) depends 

fundamentally on one or more of the system damping constants 

c , c , c or c . 

123 a 

Turning now to the high-frequency asymptotic behavior of 
Qu(oj), we have by retaining only the term proportional to oj 6 in 

D(w) and the term proportional to to 1 * In N^(io) + N 2 -j-(oj) in 

Eqs . (2.57) and (2.55) respectively as oj -*• °°, 


R(C t )" 1 (I +1 ) 

TT J.TT ^ J- 1 2 

II 21 I i I 2 t ° 2 


(2.63) 


where we have again used Eq. (2.51). From Eq. (2.58), we see that 
the asymptotic high frequency behavior of Qjj( u ) therefore can be 
expressed as 
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Q II (“) “ - 


u -v OO , 


(2.64) 


I;; 


r 2 (c 


i+- 


: x m 2 

TO -* 


Thus, as the magnitude of the transfer function Q^(oo) 

approaches unity. Moreover, for large enough oi, the deviation of 
Qjl Cm) from -1 is independent of system damping and, in our case, 

depends primarily on R 2 (C T )“Vl since I >>I . Hence 3 referring to 

2 12 

Eq. (2.60 ) , we see that as w-»-0, no dynamic forces arise in the 
test gear mesh I from the static transmission error 8-j- of that same 

mesh 3 whereas as the inertias I and I prevent all vibratory 

motion of degrees of freedom 0^ and 0 in the test rig shown in 

Fig. 1 resulting in " full transfer " of the static transmission 

error to dynamic forces in mesh I. 

Let us turn now to the transfer function jj( u ) which con- 
trols the forces generated in mesh I by the static transmission 
error in mesh II. To examine the asymptotic low frequency behavior 
of Qj j-r(w), we retain the term proportional to oi in D(oi) as given 
by Eq. 12.57) and the term proportional to w 2 in N^-^oi) + N 2 II (cd) 

as given by Eq. (2.56). According to Eqs. (2.52) and (2.59), the 
retention of these two terms gives as oi+O, 


R 2 (C tt ) ^(c +c +c ) c +c c ]gi 2 
Q t tt (u) SI 1 J 3 _ a 23 

i ( c +c +c )R 4 (C t ) -1 ( C TT ) “ 1 co 
123 1 II 


= -i 


c c 

C + 2-^ 

a c +c +c 

12 3 


01 

R 2 (Cj ) “ 1 


oi -*■ 0 


(2.65) 


Thus, as oi-»-0, Qj jj(w) also tends to zero as did Q II (oi). Moreover, 

the right-hand side of Eq. (2.65) depends fundamentally on the 

system damping constants c , c , c , and c , and is very similar 

123 a 

to the right-hand side of Eq. (2.62) which governs the behavior of 
Qjj(m) as cj-»- 0 . For c^=0, the second term in the right-hand side of 

Eq. (2.65) vanishes. 


To study the high frequency asymptotic behavior of Q 1 II (oi), 

we retain the term proportional to oi 6 in D(oi) and the term pro- 
portional to oi 4 in N 1II (oi) + N 2II (a>). According to Eqs. (2.52), 

( 2 . 56 ), (2.57), and (2.59), we then have for the asymptotic behavior 

°f Qj jj(.oi) as 
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II 

i* 


R 2 (C tt )" 1 

Q I II (0J) ~ ■ — (2.66) 

i 

which depends on the parameter R 2 (C ) _1 /I and tends to zero as 

Eenae a referring to Eq. (2.60) , we see that in both limiting 
oases as oj -»-0 and 0 )-*-“., no dynamic forces ariose in the test gear 
mesh I from the static transmission error Cjj of the slave gear 
mesh II. For intermediate frequencies 0<co<°°, however, finite 
forces will arise in mesh I from the static transmission error 
in mesh II. 


29 



THE RESPONSE 


We shall now combine the excitation Fourier series repre- 
sentation with the system transfer functions to predict response 
metrics of interest. 

Mesh and Bearing-Support Dynamic Load and Spectrum 

Since the test gears shown in Fig. 1 are assumed to move in 
(unsteady) pure rotation, the time-varying force magnitudes W^(t) 

in the test gear mesh must equal the time-varying force magnitudes 
in the bearing supports, as illustrated in Fig. 1. Recalling from 
Eq. (2.22) that x /R is the time-average force W in mesh I, we 

can express W^(t) as 

Wj(t) = x a /R + 6W I (t) , (3.1) 


where 5W^(t) represents the deviation of W^(t) from its time- 
average value x /R. From Eq. (2.60) and the linearity of the 
Fourier transformation operation, we see that 6W^(t) can be de- 
composed into two components — a component 6W-j- I (t) arising from 

the static transmission error excitation from the test gear mesh I, 
and a component 5W^ -^(t) arising from the static transmission error 

excitation from the slave gear mesh II. Hence, we may write 
Eq. (3.1) as 


w I (t) = 


V R 




+ sw 


I II 


(t) 


(3.2) 


Since in our predictions we are taking into account only the tooth 
meshing harmonics of the static transmissions error of meshes I 
and II, it follows from Eq. (2.60) that 6W^^(t) is a periodic 

function with period equal to the tooth meshing period of mesh I 
and 6Wj jj(t) is periodic with period equal to the tooth meshing 

period of mesh II. Thus, 5Wj^(.t) and 5W^ ^(t) ma y represented 

by complex Fourier series, 

°° ' i27rnf T t 

<5Wn(t) = a wjj( n ) e ( 3 - 3 ) 

n=-°° 

and 

~ i27rnf TT t 

jj(t) - ) ' jj(n)e, ( 3 • ^ ) 

n=-oo 
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where fj and f ^ are the tooth meshing fundamental frequencies of 
meshes I and II respectively, which can be expressed as 

f^ = Nj(rpm)/60 (3-5) 

and 

f^-j. = N II (rpm)/60. (3-6) 


Nj and Njj are the numbers of teeth on the test and slave gears 

respectively, and rpm denotes the number of revolutions per minute 
of all shafts in the test rig. a wil^ n ^ and a wi II^ n ^ den °te comp- 
lex Fourier series coefficients. 


It follows from the relationship between the Fourier trans- 
form and Fourier series delineated in Appendix G of Reference 4 
[especially Eq. (G6)] that we may express the complex Fourier 
series coefficients of 6W^j(t) and SW-j. ^j(t) with the aid of 

Eqs . (2.60) and (3.2) through (3.4) as 


a WII (n) = 


( c I )" 1 


Q II Ca nI )a mnI 


(3.7) 


and 


a WI II (n) ~ (C I ) lQ I Il^nll^mnll* 


(3.8) 


where 

co T — 2'rrnf T (3-9) 

nl I 

“nil = 27Tnf n> (3.10 


and where a ^ and a mn jj denote, respectively, the complex Fourier 

series coefficients a given by Eq . (1.5) for the static trans- 

mn 

mission errors arising from meshes I and II. As Eq . (1.5) in- 
dicates, a T and a is each the sum of three components that 

* mnl mnll 

are the Fourier series coefficients of the three terms on the 
right-hand sides of Eqs. (2.22) and (2.23) respectively. The 
transfer functions Q^^co) and jj(“) are defined by Eqs. (2.58) 

and (2.59) respectively. 


The temporal behavior of (t) and fiWj jj(t) can be con- 

structed from the Fourier series coefficients, Eqs. (3-7) and 
(3.8), using Eqs. (.3.3) and (3.4). These quantities are then 
added to t /R to yield the temporal behavior of the mesh and bear- 
ing support forces W^t) as indicated by Eq. (3.2). 
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Finally, assuming that the "phase” between the teeth of 
meshes I and II is unknown, we may express the one-sided power 
spectral density of W^(t) as 

oo 

$ WI (f) = (V R ) 26 ( f ) + S 2|a WII (n)| 2 5(f-nf I ) 

n=0 

OO 

+ I>Ki II (n)| 2 6(f-nf II ) , (3.11) 

n=0 

where we have been aided by Eq. (109) of Reference 4. In 
Eq. (3-11), <5 ( * ) denotes the Dirac delta function; hence, 

Eq . (3-11) describes a line spectrum with a dc contribution aris- 

ing from the mean force x /R and additional contributions arising 

from the tooth meshing harmonics nfj of the test gear mesh I and 

nf II S - Lave S ear mesh II. Notice from Eqs . (2.62) and 

(2.65) that Q II (m) and ^(.w) both vanish at co=0, so we see from 

Eqs. (3-7) and (3-8) that no dc contributions arise from the terms 
n=0 in the two sums in Eq . (. 3 . 11). 

Dynamic Force on a Single Tooth 


We turn now to obtaining an expression for the time-varying 
root stress on a typical tooth of the test gears. We first re- 
quire the time-varying force on a typical tooth. Let us denote 

this force by W^j(t). Here, W denotes force, I denotes that the 

test gears have been designated as mesh I, j denotes tooth number 

used here as a reminder that W^. ( t ) represents the force on a 

single tooth only, and the dagger denotes that here we are con- 
sidering time to be our independent variable rather than the 
lineal variable x = R0 used earlier — e.g., in Eqs. (1.2) and 
(1-3). For a fixed speed of rotation there is, of course, a 
one-to-one correspondence between x and t that is determined by 

x(t) = R0(t) = R0(t-t ), where t is a fixed reference time. 

0 0 

Let K„.(t) denote the stiffness of tooth pair j as defined 
^ 3 

by Eq. (1.16) except that here time is our independent variable. 
Stiffness Kg . (t ) varies with time because of the temporal variation 
of the line of contact on the tooth faces as the gears rotate. 


Let ut ^ (t) and uT^ ; (t) denote the elastic deformations of tooth 
3 3 

pair j of gears (1) and (2) of test gear mesh I. Then it follows 
directly from Eqs. (5) and (,8b) of Reference 4 and Eq. (I.l6) 

above that for spur gears with u^ 1 ^ and u^ 2 "* independent of axial 
location y we have 


,+ (2) 
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w ij (t) = 


K 


so 


(t)[u! (1) (t) + ut (2) (t)]. 


(3-12) 


Let denote the dynamic transmission error of test gear 

mesh I. Subscript m denotes that we are considering here only 
the effects of the mean profile modification as described in 
Reference 4 and Section 1 of the present report. The dynamic 

transmission error ? + T (t) is identical with the quantities de- 

scribed by Eqs . (2.11) and (2.13) except that here t is the in- 

dependent variable. Let m^. } (t ) and m^. ; (t) denote the mean 

profile modifications of teeth j on gears (1) and (2) of the test 
gear pair. These modifications are identical with the definition, 
Eq. (1.24), except that here time is again considered to be the 
independent variable as denoted by the dagger superscripts. Then, 
according to Eqs. (5) and (6) of Reference 4 when no random tooth 
face errors are present, we have 


ut (1) (t) + uj (2) (t) = ? m I ( t) _m sj 1) (t)-m^ 2) (t) , 

which, when substituted into Eq . (3.12) gives 


(3.13) 


W+.(t) = K+.(t) 


/ j_ \ +(1) S t(2) , . v 

S ml (t ' ) m Sj (t) m Sj (t) 


(3.14) 


Equation (3.14) is the spur gear counterpart to Eq. (8c) of Ref- 
erence 4 when no random tooth face errors are present. 


Rewriting Eq . (2.13) using our present notation gives for 

the dynamic transmission error 


4l (t) = c i W i (t) + 4 (t) ’ 


(3.15) 


where we recall that W^(t) is the total force transmitted by the 

test gear mesh I and C-J-(t) is the static transmission error of 

mesh I as defined by Eq. (2.12). If we again denote Fourier 
transforms of quantities by placing a caret over the transformed 
variable, we have for the Fourier transform of Eq. (3-15) 





(3-16) 


Multiplying both sides of Eq . (2.60) by and substituting the 

resulting expression into Eq. (3-16) gives for the Fourier trans- 
form of the dynamic transmission error of mesh I 
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3 


(3.17) 


Ci - [ 1 +q ii ( “>]^ 


I + Qj II (w)Cj I , (a / 0, 


’ml L ^11 

where (m) and Q ^(w) are defined by Eqs. (2.58) and (2.59) 
respectively, and is the Fourier transform of the static 

transmission error of mesh II as defined by Eq . (2.12). 

Equation (3.17) shows that the dynamic transmission error of 
mesh I is made up of three components, a direct component Zl equal 

to the static transmission error of mesh I, an indirect component 

whose transform is that represents the contribution to 

resulting from the system dynamic response to Cj, and finally, 

a component whose transform is jj(“)Cjj resulting from the 

system dynamic response to the static transmission error £jj of 

the slave gear mesh II. From Eqs. (2.62) and (2.65), we have 
found that Qjj( w ) and Qj jj(oa) both tend to zero as w+0 so that 

for sufficiently small to it follows from Eq. (3.17) that and 

Cj coincide — that is, the dynamic and static transmission errors 

of mesh I become identical for very low frequencies of excitation 
(as expected). On the other hand, as we have found from 

Eqs. (.2.64) and (2.66) respectively that Qjj(w) -> -1 and 

Qj jj(oj)-»- 0, so that for sufficiently large co it follows from 

Eq. (3.17) that the dynamic transmission error of mesh I tends to 

ml 

vanish. This is the mathematical manifestation of the fact that, 
for sufficiently high frequencies of excitation, the inertia of 
the gears prevents the static transmission errors of meshes I and 
II from causing any substantial lack of uniformity in the rate of 
rotation of the test gears. 

Let us now define 


it 

’mil = 


A 


1 +Qjj (10) 


( 3 . 18 ) 


and 


SnI II = Q I II (co) ^II 3 


(3.19) 


which, respectively, are the Fourier transforms of the contributions 
of the static transmission errors of meshes I and II to the Fourier 
transform of the dynamic transmission error t' T of mesh I — i.e.. 


’ml 


^ *f* ^ 4* ^ + 

? ml = ? mll + ? ml II 


( 3 . 20 ) 
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Then, denoting the inverse Fourier transforms of and 

Cl ll» res P ectivel y> by 


and 




= F _ 1 |[l+Q 


xi«»n£i} 


( 3 . 21 ) 


q,! n<« - 


1 [r T ] 

Lt ml II J 


" F ^ Q I 


(3.22) 


where we have used Eqs . (3*18) and (3.19), we then have by sub- 
stituting the inverse transform of Eq. (3-20) into Eq. (3.1*0 


W^(t) = K+J (t) [s + II (t) + C 


t 

ml II 


(t)-mj 4 (1) (t)-m^ (2) (t)] 


Sj 


Sj 

rt 


(3.23) 


Equation (3.23) expresses the dynamic force W' (t) trans- 
mitted by a typical tooth pair j of the test gear mesh I in terms 

of the tooth pair stiffness K„.(t) of the same two teeth, the 

+ 4- 

dynamic transmission error contributions ? m jj( t ) and C m j jj(t) to 

mesh I from the static transmission errors of meshes I and II re- 

t(l) + ( 2 ) 

spectively, and the profile modifications m„ : (t) and m q : (t) of 

+ (1 ) +(2) 

the two meshing teeth in contact. Quantities m^ : (t) and m„. (t) 

are considered positive when they are "equivalent" to removal of 
material from perfect involute surfaces . 

*f* 

The dynamic transmission error components C -j-j(t) and j^-(t) 

are to be evaluated using Eqs. (3-21) and (3-22). To carry out 

^ T 

this evaluation we recognize that C T contains only the tooth 

meshing harmonics of mesh I; hence, we can represent £ ^j(t) by 
the complex Fourier series 


C T-r(b) = 23 a r TT (n)e 

mil n=-°° £mll 


i2 , rrnf yt 


(3-24) 


where f^ is the tooth meshing fundamental frequency of mesh I to 

be evaluated using Eq. (3*5). Similarly, Cjj contains only tooth 

meshing harmonics of mesh II; hence, we can represent -j-j(t) 

by the Fourier series 
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? ml II (t) ^ a £ml n (n)e 
n=-°° 


i2Trnf 


II 


(3.25) 


where f^ is to be evaluated using Eq. (3-6). The complex Fourier 
series coefficients a ^ m jj( n ) and a ^ m I II^ n ^ in Eqs • (3.24) and 
(3.25) are to be evaluated from the Fourier series coefficients 

1 i 

of the static transmission errors C-j. and C j j of meshes I and II 
with the help of Eqs . (3.18) and (3-19) — i.e.. 


and 


a 


Cml II 


(n) = Q. 


II^“nII^ a mnII* 


( 3 . 26 ) 

(3-27) 


where and are defined by Eqs. (3*9) and ( 3 . 10) respec- 
tively, the a n=0 *42 , • • • are the Fourier series coefficients 

defined by Eq. (1.5) of the static transmission error of mesh I, 
and the a „ are the Fourier series coefficients defined by 
Eq. ( 1 . 5 )or the static transmission error of mesh II. Each of 
these Fourier series coefficients is the sum of three terms as 
indicated by Eq. (1.5). 

The left-hand sides of Eqs. (3-26) and (3-27) are the com- 
plex harmonic amplitudes of the contributions to the dynamic 
transmission error of mesh I from the static transmission error 
harmonic amplitudes of meshes I and II respectively. The right- 
hand side of Eq. (3-26) exhibits a "direct self" contribution 
a _ from the static transmission error plus the dynamic incre- 

ment contribution Qjj( w n j) a mn i > whereas the right-hand side of 

Eq. (3-27) Qj 11 (“nil ^ a mnII the i ncrement to mash I from 

mesh II which arises from dynamic effects only. (We have shown 
earlier by Eqs. (2.62) and (2.65) that Q^^ui) and ^(w) both 

tend to zero as oo->0. Hence, for i = I or II and j = I or II, 

Q..(w) represents the dynamic contribution to the (dynamic) 
a J 

transmission error complex harmonic amplitude at mesh i from a 

unit amplitude harmonic component of the static transmission 

error at mesh j . This fact suggests that the term "transmission 

error dynamic increment function" or the more brief term " dynamic 

increment function" be applied to the transfer functions Q..(oj). 

aj 

This concept, and the foregoing analysis, readily generalizes to 
gearing systems with an arbitrary number of dynamically interact- 
ing meshes. 
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Equations (3.23) through (3*27) are one set of final expres- 
sions for evaluation of the temporal behavior of the forces 

^lj(t) transm itted by a typical pair of teeth on the test gears. 

Decomposition of dynamic single tooth force into static and 
dynamic increment contributions . Instead of using the form of 

4 - 

W T .(t) given by Eq. (3.23), we can again begin with Eq. (3.1*0 
J -f- 

and decompose the dynamic transmission error ^ -j-(t) of mesh I into 

the original static transmission error £ ^(t) of mesh I plus 

dynamic increment contributions to mesh I arising from meshes I 
and II — I . e . j 

5 ml^ = ? ml (t) + 6S mII (t) + 6C mI II (t) s (3*28) 

where £ j(t) is the static transmission error of mesh I defined by 

Eqs. (1.3) through (1.5), and and <5?^j jj(t) are the dynamic 

increments to the transmission error of mesh I arising from 
meshes I and II respectively. Then, by substituting Eq. (3-28) 

•j* 

into Eq. (3-14), we see that the dynamic force W T .(t) can be 

.j- 

decomposed into components W^j o (t) and 6W^^. (t) arising, respectively 
from the static transmission error £ m j(t) and the dynamic incre- 
ments <5£^n(t) + 6^1 ij(t) - i.e., 

W|j(t) = W+ (t) + 6W+ (t) , (3.29) 

where 

Wjj Q (t) ^ 4j(t)[? m l(t)-4j 1)( t)- m Sj 2) (t)], (3-30) 

and 

( t ) £ K Sj (t)[6 4ll (t)+6 4l II (t)] ‘ (3-31) 

The static transmission error ? m j(t) can be expressed in terms of 

its Fourier series coefficients a T by 

mnl 

°° i 27 rnf T t 

VE (t) = E a mnl e * (3.32) 

n=-oo 
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where f_ is the tooth meshing fundamental frequency of mesh I to 
be evaluated using Eq . (3-5), and the a mn j are the Fourier series 

coefficients of the static transmission error of mesh I defined 

by Eq. (1.5). The dynamic increment 6c^jj(t) can be represented 

in terms of its Fourier series coefficients Qj j ( w n j ) a mn j by 


5? 


t 

mil 


(t) 


E Q II<“nI 
n=-°° 


)a T e 
mnl 


127 rnf jt 


(3.33) 


and the other dynamic increment 6 j- jj(t) can be represented by 




t 

ml II 


(t) = £ Q 

n=-°° 


i2irnf-j-.j.t 
tt (oj TT )a T T e, 

I II nil mnll s 


(3-34) 


where and ^ nI j are defined by Eqs. (3-9) and (3-10) respec- 

tively, f^ is the tooth meshing fundamental frequency of mesh II 
to be evaluated using Eq. (3-6), and the a mn ]-j are the Fourier 

series coefficients of the static transmission error of mesh II 
defined by Eq. (1.5). The transmission error components in Eqs. 
( 3 . 23 ) and (3.28) are related by 

+ • <3.35) 


and 


.+ 

’ml II 


(t) s Sc 


t 

ml II 


(t) 


(3-36) 


The representat ion of the dynamic force on a single tooth 
by Eqs. (3*29) through (3-34) has the advantage of separating out 

the force increment SW-j-j(t) arising from the dynamic response of 

the gearing system. Notice from Eqs. (3-33) and (3-34) that 

QllCoini) and Q T TT (m inTT ) govern the dynamic response contributions 


‘I II s nil 


from the static transmission errors of meshes I and II to the 
force increment 6W^.(t). Furthermore, since Qjj(0) = 0 and 


Q 


I II 


(0) = 0, we see from Eqs. (3.31), (3.33), and (3-34) that 


when all a mn j and are zero except for n = 0, we have 

SWjj(t) = 0. Eence , when the fluctuating components of the static 
transmission errors of meshes I and II are zero, Eqs. (3.31), 
(3.33), and (3.34) imply that no force increment Swt.(t) on a 

J 
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f 


single tooth arises from the dynamic response of the gearing 
system. 

Finally, we note that a somewhat more accurate expression 
for the total mesh force W^(t) than that given by Eqs . (3-2) 

through (3-10) can probably be obtained by summing the indivi- 

4 * 

dual tooth forces W-j-j(t) given by Eq. (3-29) over all teeth in 

contact at any given instant of time. However, in carrying out 
this summation, care must be exercised to take into account the 
appropriate "time delays" of the individual tooth forces. 

Root Stress on a Single Tooth 

Cornell [8] has modified a formula developed by Heywood [9] 
for the root fillet tensile stress in a gear tooth. Additional 
comments by Heywood and others on this problem appear in the 
discussion of the paper by Kelley and Pederson [10]. Cornell’s 
modification of the Heywood formula is readily adapted to the 
geometry of gear teeth of nominal involute design — his modifica- 
tion can be expressed as 



where a is the root stress and W. is the tooth loading denoted by 

J 

Wjj in the case of mesh I, and given by Eqs. (3-23) or (3-29) - 

The remainder of the notation in Eq. (3-37), and expressions for 
the parameters contained therein, are contained in Appendix B of 
this report. The procedure developed by Cornell [8] for determina 
tion of the position of the maximum root fillet tensile stress 
also is described in Appendix B. 
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PROFILE DESIGN FOR MINIMUM EXCITATION 


Profile Modification to Compensate for Tooth and 
Gear Body Deformations at the Design Loading 


Let c(x) denote the (dynamic) transmission error as defined 
in Reference 4 and expressed by Eq. (2.6) of this report. Here, 
the transmission error is expressed as a function of the coordi- 
nate x, which can be related to the angular position 0 of either 
gear of a meshing pair by Eq. (1.2). The requirement that the 
teeth of the gear pair remain in contact at each value of x 
dictates [4] that the tooth surfaces must move together relative 
to the tooth surfaces of rigid perfect involute gears by the 
distance 


4 ( x ) = u^Cx) + u? 2) (x) + m^(x) + m? 2) (x) , 

J J J J 


(4.1) 


where u[ (x) and u^ (x ) are the elastic deformations of the jth 

2 2 ('q') (?) 

(meshing) teeth on gears (1) and (2), and m: '(x) and m: ; (x) are 

2 2 

the modifications of the profiles of the same teeth on gears (1) 

and (2) from perfect involute profiles. Quantities £, u^*^ and 
( * ) 

m are measured in a direction defined by the intersection of 

the plane of contact and the transverse plane. Quantities u[*^(x) 
( . ) J 

and m. (x) are defined as positive when they are "equivalent" to 
J 

removal of material from the surfaces of unloaded 3 perfect j uni - 
formly spaaed involute teeth [4]. Thus, a positive value of £(x) 
corresponds to tooth surfaces coming together relative to those 
of rigid perfect involute gears. 


It follows from the definition of the static transmission 
error that if a real pair of meshing gears is to transmit exactly 
uniform angular motion then £(x) must be a constant [11-13], In 
addition, however, if we wish to insure that the tooth faces under 
loading take on the form of perfect uniformly spaced involute teeth, 
we must require that the tooth faces from each gear be modified to 
compensate for their own elastic deformations — i.e., that the geo- 
metrical deviation of each tooth face take on the form 

m j 1 ^ ( x ) = —u ^ 1 ^ ( x ) (4.2a) 

*J U 

m (2) (x) = -u? 2) (x) (4.2b) 

J J 

which will yield zero static transmission error, ?(x) = 0, when 
combined with Eq. (4.1). When Eqs . (4.2a) and (4.2b) are both 
satisfied, we are guaranteed that all tooth pair contact remains 
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II 

if 


in the plane of contact; however, satisfaction of the single 
condition 


(x) 

J 


m< 2) (x) = 

J 


uS 1 * ( x ) +u ^ 2 ^ (x) 
3 3 


(4.3) 


will not guarantee this behavior, and therefore, will not guaran- 
tee the transmission of exactly uniform angular motion [l4j. 

To evaluate the right-hand sides of Eqs . (4.2a) and (4.2b), 

we require the loading W.(x) carried by each tooth pair j in 

J 

contact, as well as the appropriate stiffness characterizations. 
The tooth pair loadings W.(x) must be constrained by the relation- 
ship : 4 


W = £ W, (x) 

d j J 


(4.4) 


where is the total loading transmitted by the gear pair, where 

the summation over j in the right-hand side includes all tooth 
pairs in contact, and where all forces in Eq. (4.4) are in the 
direction defined by the intersection of the plane of contact 
and the transverse plane. The constant force is the design 

loading of the gear pair for minimum vibration excitation. Follow- 
ing Harris [ll, 12] and Remmers [13], our design approach will be 

to specify the individual forcing functions W.(x) so that the 

J 

overall loading constraint, Eq. (4.4), is satisfied. In particular, 
this design approach permits us to control the individual tooth 

loading functions W.(x) at the initiation of tooth contact. 

J 

( • ) r t 

Let c:.;(x) denote the influence function L 1 5 J that describes 
0 J 

the surface deformation of the contact patch of tooth j on gear 
(•) due to a "unit" loading of the contact patch of tooth j* on 
the same gear, as illustrated in Fig. 2. The center positions of 

( • ) 

the contact patches are determined by the rotational position 0 

of gear (•) as described by the lineal variable x = R^ ^ 6 ^ ^ , 

where all superscripts (•) refer to the same gear of the pair. 

The direction of all such deformations and loadings is defined by 
the intersection of the plane of contact and the transverse plane. 
The width 2b of the contact patch on tooth j and the semi-elliptical 
pressure distribution under this contact patch are determined by 

the prescribed loading W.(x) transmitted by tooth pair j and the 

J 

two local radii of curvature of the jth teeth on gears (1) and 
(2) at their contact point — as determined, e.g., from Eqs. (235) 
and (23 6) on p. 4l8 of Reference [l6]. Thus, the influence func- 
tions c^.)(x) include bending, shear, Hertzian contact, and gear 
J J 
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body deformations, where these latter deformations are particularly 
important in determining the deformation of one tooth caused by 
loading on another tooth, which is the case j ’ f j illustrated in 

Fig. 2. The influence functions c..,(x) can be computed using 


J J 


finite element techniques; however, particular care must be exer- 
cised to insure accurate computation of the contact deformation 
component . 


From the influence functions and the prescribed loading func- 
tions, we can compute the deformation of tooth j on gear (•) as a 
function of x by superposition 


u 


(•) 

j 


Cx) 


I 

j’ 



Cx) 


Wj , (x) 


(4.5) 


where the summation over j ’ includes all teeth in contact at posi- 
tion x = ^ 0 f ^ By combining Eq. (.4.5) written for gears (1) 

and (2) of the meshing pair with Eqs . (4.2a) and (4.2b), we obtain 
our final profile modification formulas for spur gears: 


m$ 1) (x) = - E c^J Cx) W. , (x) 

J j I J J J 

m( 2 \x) = -E c(?|(x) W . , ( x ) , 

J -? T J J J 


(4.6a) 

(4.6b) 


where the loads W.(x) transmitted by the individual pairs of teeth 

J 

are subject to the constraint, Eq . (4 4) , and where a positive 
( • ) 

value of m: (x) corresponds to " removal " of material from a perfect 

J 

involute profile. From Eqs. (4.6a), (4.6b), and (4.4), we see that 
it is possible to modify the profiles of both gears of a meshing 
pair to completely compensate for tooth deformations only at a 
single design loading W^. 


Thus far, we have ignored the thickness of the layer of lubri- 
cant between pairs of meshing teeth. These effects could be in- 
cluded, if known, by adding one-half of the local layer thickness 
to the right-hand sides of Eqs. (4.6a) and (4.6b). 


Discussion : When the profile modifications 


m^ 1 } ( x ) 

J 


and 


m[ 2 ) (x) 
J 


of all teeth on a pair of meshing gears are machined to exactly 
satisfy Eqs. (4.6a) and (,4.6b) , the static transmission error of 
that gear pair is set equal to zero when the gears are operating 
at the design loading W , . For the present discussion, let us assume 
that the profile modifications of the test gears, shown as mesh I 


43 



in Pig. 1, have been machined so as to satisfy Eqs. (4.6a) and 
(4.6b). In this case, all Fourier series coefficients 
a -r n = ±1 , ±2 , • • • of the static transmission's error are identi- 

cal to zero. It then follows from Eqs. (3*3) and (3-7) that the 
time varying component 6W II (t) of the mesh and bearing force 

attributable to mesh I in Eq. (3-2) is identically zero. Similarly, 
if the teeth of the slave gears in mesh II were modified according 
to Eqs. (4.6a) and (4.6b), the time varying force component 
SWj jj(t) in Eq. (3-2) attributable to the slave gears would be 

identically zero as can be seen from Eqs. (3*4) and (3*8). 


Turning now to the dynamic force components on a single tooth, 

•j- 

we recall that 6W^j(t) given by Eq. (3*31) is the component of the 

force history on a single tooth of mesh I attributable to the 

vibratory excitation from meshes I and II. Examination of 

Eqs. (3-31), (3.33) and (3-34) shows that when the teeth of the 

gears in meshes I and II are designed according to Eqs. (4.6a) and 

(4.6b) to eliminate the static transmission errors of those meshes, 

their corresponding contributions to 6wt.(t) vanish thereby eliminat- 
f ij 

ing <5Wjj(t). Furthermore, we see from Eq . (3.32) that ? m j(t) also 
vanishes in this case, thus leaving 


w L. (t) = * K s 3 (t) 


t (1) / , x , +(2) f N. 

m s j ( t ) + m S j (t) 


(4.7) 


as the time-varying force history on a typical tooth of mesh I. 


Each of the single tooth loading components in Eqs. (3.29) 
through (3-31) gives rise to a corresponding component of root 
stress. Thus, when the static transmission errors of the mesh 
and slave gears are eliminated at a design loading W^, the com- 


ponents of root stress arising from the single tooth loading 

component 6W^. (t) of Eq. (3-31) are eliminated, and the resultant 

root stress is that resulting from the tooth loading component 
given by Eq. (4.7). 


Families of Tooth Loading Functions 


The design procedure suggested in the previous section is to 
modify the profiles on the teeth of two meshing gears as pre- 
scribed by Eqs. (4.6a) and (4.6b), where the loading histories 
Wj(x) on the individual teeth must satisfy the constraint given 

by Eq. (4.4). This procedure is most easily carried out by speci- 
fying the loading histories W (x) subject to the constraint, Eq . (4.4), 
and then solving for the profile modifications given by Eqs. (4.6a) 
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and (4.6b). We recall according to Eq. (1.2) that independent 
variable x is the roll distance 


x = r^ 1 ^ 1 ^ = R^e*' 2 ) 


(4.8) 


( • ) ( • ) 

which is common to both meshing teeth, and R fe and 0 are the 

base circle radius and the roll angle of each gear of the meshing 
pair where (•) designates the gear (1) or (2) of the pair as 
appropriate . 

Initially, let us restrict our considerations to spur gear 
pairs with contact ratios between one and two. There exists an 

infinite number of individual tooth loading histories W.(x) that 

J 

satisfy Eq. (4.4). Figure 3 illustrates the individual tooth load- 
ing histories W.(x) and the manner in which they add up to satisfy 
J 

the constraint, Eq . (4.4). The ordinate of Fig. 3 is force trans- 
mitted between the teeth of driving and driven gear in the direction 
defined by the intersection of the plane of contact and the trans- 
verse plane. The abscissa in Fig. 3 is roll distance x = R^0- 

Alternating intervals of single tooth pair contact and two tooth 
pair contact are shown in the figure. The, heavy line represents 
the loading and unloading force history W. (x) on the jth tooth of 
one of the pair of meshing gears. The dashed lines Illustrate the 
force histories on teeth j-1 and j+1 as indicated. The roll dis- 
tance interval during contact of a particular tooth is L and the 
base pitch is A. 

In order for the condition indicated by Eq . (4.4) to be satis- 
fied during regions of single pair contact, the design loading 

must be carried by that tooth pair as indicated in Fig. 3. On the 
other hand, during regions of two tooth pair contact, the total 
loading W^ must be shared by the two tooth pairs in contact. This 

requirement means that distances A and A' shown in Fig. 3 must be 
identical . 

If the direction of rotation of the gear pair were to be re- 
versed, and if in this reversed situation the former driving gear 
becomes the new driven gear and vice versa, the "loading side" of 
the teeth in the former situation also remains loaded in this re- 
versed situation. This reversed situation in effect reverses 
the direction of the x-axis in Fig. 3. This reversal of the direc- 
tion of the x-axis suggests that it is reasonable to require the 
tooth loading histories W.(x) to be symmetric about the midpoint 

x of the tooth contact interval L — i.e., to require that 
o 

W.(x-x ) = W.[-(x-x )]. (4.9) 

J 0 . ■ J 0 
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FIG. 3. TOOTH LOADING FUNCTION FOR CONTACT RATIO BETWEEN 
ONE AND TWO. 
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II 


In particular, this condition implies that the abruptness of loading 
onset and loading release on a given tooth be the same. Whenever 
the condition of Eq. (4.9) is met, one can readily see from the 
periodicity of the Individual loading curves in Pig. 3 that dis- 
tances A* and B must the same; hence, the symmetry requirement 
Eq. (4.9) together with the loading constraint Eq. (4.4) jointly 
Imply that distances A and B In Fig. 3 be the same. In other words , 
the single tooth loading curve in a region of two tooth pair contact 
must be an odd function about the axis W = W^/2 when Eqs. (4.4) and 

(4.9) are satisfied. Finally, within each region of two tooth pair 
contact, we should like to require that practical tooth loading 
curves be monotonic , which implies that the slope of the tooth 
loading curve not change sign within that region. Sketches of 
tooth loading functions have appeared in References 11 and 13. 

The non-constant portions of tooth loading functions satisfying 
the mirror image property, Eq . (4.9), are defined entirely by their 
behavior within a single two tooth pair contact region as is easily 
seen from Fig. 3* Let us define a new independent, dimensionless 
roll distance variable 


x-x 

A Q 

- Tl=ST/2 


(4.10) 


where x^ is the center of the two tooth pair contact region as 

illustrated in Fig. 3- Thus, the new normalized roll distance 
variable ranges from -1 ^ x ^ 1 within the two tooth pair contact 
region. Let us also define a dimensionless tooth loading variable 


W. 

J 


* v ( y 2) 

* V 2 


O.ll) 


which also ranges from -1 < W . < 1 when W^. ranges from 0 to . 

A simple family of tooth loading functions that satisfy all of 
the above mentioned criteria is the family 


W ( x ) = x + 


M 


sin(Trx) , 


-1 >< x x 1 


(4.12) 


where for W.(x) to remain monotonic, we must have 0 ^ s ^ 2. The 

J 

parameter s is the slope of Wj(x) at x = -1, i.e.. 


dW. 
S = i- 


dx 


(4.13) 


x=-l . 
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Figure 4 illustrates the behavior of Eq. (4.12) over the range of 
values of s for which W.(x) is monotonic. 

For tooth pairs with profile modifications designed using 
Eqs . (4.6a) and (4.6b), the sensitivity of the actual roll distance 
contact range L to variations of the actual mesh loading from the 
design loading depends strongly on the slope of profile modifica 

tion at the points of tooth contact initiation and termination. 

This profile modification slope, in turn, depends strongly on the 
slope of the tooth loading functions W.(x) at points of contact 
initiation and termination which is governed by the dimensionless 
slope s defined by Eq. (4.13). [Recall that s also controls the 
slope at contact termination because of the mirror image property 
Eq. (4.9)]. Thus, if it is desirable to minimize changes in the 
actual value of L over a wide range of mesh loadings, then a large 
value of s should__be_used in the profile modification design 
loading function ¥.(x). We can see from Eq. (4.12) or Fig. 4 that 

the sinusoidal-based loading function is limited to initial dimen- 
sionless slopes s less than or equal to 2 if it is to remain mono- 
tonic. 


A second family of tooth loading functions that does not have 
this limitation is the polynomial-based normalized family 

Wj (x) = x + * V > -1 « x < 1 (4.14) 

where the parameter r is permitted to take_on_only the odd integer 

values r = 3,5,7,' ,# . This definition of W.(x) remains monotonic 

1 

provided that the initial (dimensionless) slope £ remains within 
the range 0 << s ^ r. Thus, if a large value of s is desired, a 
correspondingly large value of r must be chosen in the loading 
function, Eq . (4.14). 

Figures 5 and 6 illustrate the normalized families of tooth 
loading functions of Eq . (4.14) for values of r = 3 and 7 respec- 
tively. From these two figures, we see that as r is increased 
the loading function of Eq . (4.14) permits an ever increasing 



character. Figure 7 shows h£w the loading functions of Eq . _C_4 .14) 
vary when the initial slope s is held to the constant value s = 3, 
while r is varied from values of 3 to 9- 

Figures 4 through 7 collectively illustrate the considerable 
flexibility of shapes available in the two families of loading 
functions described by Eqs. (4.12) and ( 4 . 1 4 ) . Notice that 
Eqs^ (4. 1_2 ) and_(4.l4) both reduce to the same linear function 
W.(x) = x when s is taken to be unity. 

J 
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FIG. 4. NORMALIZED S I N USO I DAL -BASED FAMILY OF TOOTH 
LOADING FUNCTIONS. 






NORMALIZED TOOTH LOADING, W 





NORMALIZED TOOTH LOADING, W = 



NORMALIZED ROLL DISTANCE, x 


FIG. 7. NORMALIZED POLYNOMIAL-BASED FAMILY OF TOOTH 
LOADING FUNCTIONS FOR 7=3. 
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Families of tooth loading functions satisfying Eq. (4.4) for 
contact ratios larger than 2 can readily be generated. Figure 8 

illustrates a linear tooth loading function W.(x) for a contact 

J 

ratio between 2 and 3 that satisfies Eq. (4.4). Here a in addition 
to the parameters L and A, and the mesh loading W^, we have an 

additional free parameter F, where FW^ is the maximum value of 

W.(x), as can be seen from Fig. 8. For the portion of the tooth 
J 

loading curve W.(x) to remain monotonic in the region x x , the 
J o 

parameter F must be chosen to lie within the range 1/2 F ^ 1 . 

In practice, F can be chosen to be the value that minimizes the 

maximum root stress of a tooth. 

Examination of each small, dashed, rectangular box in Fig. 8 
shows that the sum of the two individual tooth loading curves within 
each box equals the height of the box. However, this is exactly the 
requirement placed on the tooth loading curves within each two-tooth 
contact region illustrated in Fig. 3- Thus, each straight line seg- 
ment and its mirror image in Fig. 8 could be replaced by any of the 
loading curves shown in Figs. 4 through 7 and their respective mirror 
images to yield various families of curved tooth loading functions — 
all of which would satisfy the fundamental design constraint 
Eq. (4.4) and the symmetry property, Eq . (4.9). In particular, use 
of the polynomial-base family illustrated in Figs. 5 through 7 
permits, for any choice of the parameter F shown in Fig. 8, a com- 
pletely free choice of the slope of W.(x) at contact Initiation 

J 

and termination provided that r is chosen sufficiently large. For 
a given choice of contact initiation slope, minimization of the 
maximum root stress on a tooth might require simultaneous varia- 
tions of the two parameters r and F in the stress minimization pro- 
cedure. (.Hertzian) contact stress minimization might also be taken 
into account in such a minimization procedure. 

Excitation Spectra of Tooth Loading Functions 

A tooth loading function, such as that illustrated by the 
heavy line in Fig. 3, has the potential for exciting vibratory 
motion of the tooth itself [171. Such potential tooth vibrations 
are more likely to be a problem in high speed gearing, since an 
increase in gear rotational speed will reduce the total duration 
during which a tooth is loaded, thereby "shifting" the spectrum 
of the tooth loading function to the higher frequency range where 
tooth resonance frequencies are found. To assess the potential 
for such excitation of tooth resonant frequencies, the Fourier 
transforms of the tooth loading functions W.(x) are required. 

Each of the dimensionless loading functions'] Eqs. (4.12) and (4.14), 
can be combined with its mirror image to determine a dimensional 
tooth loading function as illustrated in Fig. 3- The Fourier 
transform of the resulting sinusoidal-based loading function with 
origin at x - x is 
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FIG. 8. 


TOOTH LOADING FUNCTION FOR CONTACT RATIO BETWEEN 
2 AND 3. 
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F [W (x);g] A [ W (x)e 12Trgx dx 

X d Jo 


= V 


sin (irAg ) 

| sin[ir (L-A )g] 

1-s (L-A ) 2 g 2 

tt Ag 

j TT (L— A ) g 

1- (L-A ) 2 g 2 _ 


(4.15a) 

(4.15b) 


whereas, the .Fourier transform of the resulting polynomial-based 
loading function with origin at x '= x n is 

x “ j ' ' 'd“ trAg \ a tt (L-A)g 


r ru / rVr i _ H a sin(TTAg) I - sin[7T(L-A) g ; 
F x LW 1 (x) ’ gJ " W d A ttAs | 3 TT ( L-A ) g 


^ l 2 


r : 


Tr^kT 


cos [tt ( L-A )g+7j-kTr] 
[TT(L-A)g] k 


r = 3,5,7, , s^r. 


(4.16) 


The right-hand sides of Eqs. (4.15b) and (4.16) both have the 
form of ¥ d A sin(TrAg)/(TTAg) times a quantity within curly braces. 

Thus, we can study the effects of the dimensionless slope s at 
contact initiation and termination [Eq. (4.13)] by examination of 
the behavior of the quantities within the curly braces in 
Eqs. (4.15b) and (4.16). The magnitudes of these quantities are 
plotted on log-log coordinates for the sinusoidal-based loading 
function in Figs. 9a to 9d , and for the polynomial-based loading 
function in Figs. 10a to lOd and 11a to "lid. The ordinate in 
Figs. 9a through lid is the (dimensionless) magnitude of the quan- 
tities within the curly braces in Eqs. (4.15b) and (4.16), and 
the abscissa is measured in units of 

q = (L-A)g = (£ -l) Ag, (4.17) 

where g is the transform variable indicated by Eq. (4.15a). 


To interpret the dimensionless frequency q, we note that since 
A is the base pitch, its reciprocal 1/A is the tooth meshing fun- 
damental frequency in the frequency units g, which we denote by 

g, , , , . Thus 

& tooth- mesh 


Ag = — S_ = S = p , (4.18) 

A" & tooth mesh 

where p is the tooth meshing harmonic number which is unity at the 
tooth meshing fundamental frequency. Combining Eqs. (4.17) and 
(4.18) gives 
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q=(L-4)g 
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FIG. 10. SPECTRUM FACTOR WITHIN CURLY BRACES IN EQ. (4.16) 
ON LOG-LOG COORDINATES FOR POLYNOMIAL-BASED TOOTH 
LOADING FUNCTION FOR r= 3. 
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SPECTRUM FACTOR WITHIN CURLY BRACES IN EQ. (4.16) 
ON LOG-LOG COORDINATES FOR POLYNOMIAL-BASED TOOTH 
LOADING FUNCTION FOR r=7. 








q 


(4.19) 



g 

®tooth mesh 



> 


where we recall that (L/A) is the transverse contact ratio which is 
generally in the approximate neighborhood of (L/A) ~ 1.5 

Examination of Eqs. (4.15b) and (4.16) shows that both equa- 
tions have the same high frequency asymptotic behavior : 


F x [Wj (x) ; g] 


A sin ( it Ag) (- sin[ir(L-A)g] | 
d A irAg | s it ( L— A ) g ) 

a sin(TrAg) i - sin(Trq) ) 
d n ' 7T Ag j irq j ‘ 


(4.20) 


¥e see that the quantity within the curly braces in Eq_^ (4.20) is 
independent of r, from which it follows that provided s / 0 the 
high frequency asymptotic behavior of all curves with the same value 
of s in Figs. 9 through 11 should be the same. The figures are in 
agreement with this behavior. Furthermore, from Eqs. (4.10), (4.11), 
and (4.13), we see that 


(4.21) 

contact initiation 


s = 


dW. 

J 


dx 


- (L-A) 

W, 


dW, 


dx 


x = -l 


where the slope of the tooth loading function is evaluated at the 
point of contact initiation. Combining this result with the first 
form of Eq. (4.20) gives for the high-frequency asymptotic behavior 
of F [W (x) ; g] : 

A J 


F x [ w j (*) ;g] 


sin(frAg) sin[Tr (L-A)g] 


2 2 
ir g 


dx 


(4.22) 

contact initiation 


whose envelope Is ( T.V )'* (dW^/dx) contact lnltlatlon - ‘Thus, we 

see that the high frequency content of the Fourier transform of 
the tooth loading function is directly proportional to its slope at 
the point of contact initiation. 


As a consequence of this fact, we have a trade-off in choosing 

the tooth loading functions W.(x) for high-speed gearing. On the 

J 

one hand, a large value of the slope of W.(x) at points of contact 

J 

initiation and termination will minimize variations in the tooth 

contact region for loadings other than the design loading. On the 

other hand, such large values of the slope of W.(x) will also give 

J 
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rise to more energy in the high frequency portion of the spectrum 

of W.(x), and therefore will be more likely to excite tooth reso- 

J 

nant vibrations. For tooth loading functions with the same slope 
magnitudes at points of contact initiation and termination, 

Eq. (4.22) provides a simple formula for estimating the high fre- 
quency region of the Fourier transform of (x) for investigating 
the importance of the resonant response of gear teeth. 



APPLICATION 


The theoretical results of the preceding sections of the 
report are applied to the NASA Lewis four-square test apparatus 
in this section. 

Profile Modification for Minimum Vibration Excitation 

The initial step in the profile modification design pro- 
cedure is to specify the individual tooth loading function W. (x) 

tJ 

illustrated in Fig. 3, where W.(x) is subject to the constraint 

given by Eq . (4.4). The simplest loading function to satisfy 

Eq. (4.4) is that characterized by the dimensionless function 

W. (x) given either by Eq. (4.12) or Eq . (4.14) for the case where 

s = 1. [For this value of s, Eq. (4.12) and (4.14) both describe 

the dimensionless linear loading characterization 

W. (x) = x, - 1 << x ^ 1.] For this case, the full dimensional 

loading function W.(x) becomes the function with linear load 

sharing curves shown in Fig. 12, which is symmetric about the 
pitch point. The abscissa of Fig. 12 represents "roll distance" 
6x = where 60 is an increment of gear rotational angle in 

radians and is the base circle radius of the same gear. The 

maximum potential roll distance over which a tooth can be loaded 
can be determined by Eqs . (D8), (Dll), and (D12) of Ref. 4 from 
the addenda radii and pitch radii of the two meshing gears and 
the pressure angle. Both test gears in the NASA Lewis four- 
square apparatus have R^ = 4.7625 cm (1.875 in) and R = 4.445 cm 

(1.750 in) for their addenda radii and pitch radii respectively. 
Their pressure angle is 20°. Utilizing Eqs. (D8), (Dll), and 
(D12 ) of Ref. 4, we find the maximum potential value of roll 
distance for these test gears to be L nom ^ nal = 1.53530 cm 

(.604450 in) as shown in Fig. 12. However, the tip rounding 
radius on the test gears reduces the actual maximum value of roll 
distance to a value of about L = 1.37450 cm (.5411406 in). Fur- 
thermore, it was our desire to allow for at least a 20% increase 
in loading of the test gears beyond the design loading, and to 
allow for tooth spacing errors without encountering tip inter- 
ference. Thus, we have chosen the design value of roll dis- 
tance to be L = 1.22971 cm (.48414 in). This value can be com- 
puted from the values shown in Fig. 12 — i.e., 

1.53530 - (2x. 152794) = 1.22971 cm. The values of roll distance 
increment during load sharing of .29235 cm are computed from the 
design value of L = 1.22971 cm and the base pitch A. The base 
pitch is computed from the pitch radius R = 4.445 cm, number of 
teeth N = 28, and pressure angle <f> = 20° by the relation 
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FIG. 12. TOOTH LOADING FUNCTION FOR COMf 
STIFFNESS [LENGTH DIMENSIONS Af 



ATION OF TOOTH 
IN CM (IN).]. 




A 


.93730 cm ( . 36902 in) . 


(5.1) 


2irR 

N 


COS<t> = 


The roll distance increment during load sharing is L-A (see 
Fig. 3.)» which in the present case is 


L-A = 1.22971 - .93730 = .29241 cm (.1151 in) 


(5.2) 


as is shown in Fig. 12. The tooth loading function shown in Fig. 12 
was constructed in the above manner. These values result in a 
design contact ratio of (L/A) = ( 1 . 22971/ . 9 37 30 ) = 1.312, which 
is a completely acceptable value. 

The tooth loading diagram shown in Fig. 12 was transmitted to 
NASA Lewis and used there to compute the tooth deformations re- 
quired to carry out the profile modification designs. See Fig. 4 
of Ref. 18, where the roll angle increments 60 in degrees are 
obtained from the roll distances 6x shown in Fig. 12 by the relation 

66 deg 2irR b 6x , (5-3) 


where, for the test gears in question, the base circle radius is 
R b = 4.17693 cm (1.64446 in). 

Reference 18 computes the deformation under the loaded tooth, but 
not the component of deformation under neighboring teeth due to 
that loading. Thus, to illustrate the profile design method de- 
scribed in the preceding section, we shall assume here that when 
only one tooth is loaded, the deformations of all neighboring teeth 
are zero. That is, for the present approximate profile modifica- 
tion design, we shall assume that c..,(x) = 0 when j ^ j ' , where 

J J 

the c..,(x) are the influence functions in the profile design 

formulas given by Eqs. (4.6a) and (4.6b). Employing this assump- 
tion, and recognizing that, in the present application, both mesh- 
ing test gears are identical, we see that both design equations 
(4.6a) and (4.6b) are reduced in this (approximate) application to 

m.(x) = -c..(x) W.(x) , (5.4) 

J JJ J 

which is the negative of the tooth deformation due to application 
of the loading W.(x) illustrated in Fig. 12. However, according 

to our sign convention, a negative value of m. (x) is equivalent to 

J 

an addition of material to the involute design. Thus, the approxi- 
mate result of Eq. (5.4) states that an amount of material depth 
equal to the tooth deformation under the loading W.(x) must be 

J 

added to the perfect involute tooth surface to minimize the vibra- 
tion excitation. This depth is measured in a direction defined by 
the intersection of the base plane and a plane normal to the gear 
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axes. This is the direction normal to the perfect involute 
surfaces before modification. 

The tooth deflections caused by the point loading function 
shown in Fig. 12 have been tabulated in Table III of Ref. 18, 
which has been reproduced as Table 1 of this report. Thus, the 
deflections listed in Table 1 are the depths of addition of ma- 
terial to the test gear pure involute profiles as described by 
the approximate result of Eq . (5.4). 

Normally, profile modifications are thought of in terms of 
removal of material from perfect involute profiles. The amounts 
of material "removal" equivalent to the material "additions" 
listed as "Deflections" in Table 1 can be obtained by subtracting 
the deflection entries in Table 1 from the maximum deflection 
entry. The column entitled "Profile relief" (in) in Table 2 was 
computed in this manner by subtracting the corresponding "Deflec- 
tions" (in) in Table 1 from 6.1733 E-4 . The values of profile 
relief in mm were then computed by multiplying the (in) entries 
by 25.4. 

Examination of the profile relief values given in Table 2 
shows maximum values of profile relief at the two ends of the 
contact region (near tooth tip and root), with relief values 
monotonically decreasing to a value of zero in the center portion 
of the tooth. This general behavior is consistent with current 
industry practice, and is a consequence of the monotonically in- 
creasing and decreasing behavior of the tooth loading function 
shown in Fig. 12. Any of the sinusoidal or polynomial-based load- 
ing functions discussed in the preceding section would yield pro- 
file relief functions with the same general behavior as that 
tabulated in Table 2. 

The teeth of both of the nominally identical test gears must 
be modified by the amounts indicated in Table 2 to compensate for 
their deflections under loading. The smaller values of roll angle 
designate the tooth regions nearer the tooth root, and the larger 
values of roll angle designate the tooth regions nearer the tip. 
The profile modification for approximate minimization of the 
vibratory excitation has the same shape as the tooth deflection 
function shown in Fig. 13. The modifications of the root and 
tip regions beyond the design zone of contact are to be obtained 
by linear extrapolation of the profile relief values shown in 
Table 2. This procedure is equivalent to a linear extrapolation 
o/f the two ends of the curve shown in Fig. 13. This increased 
relief in the tip and root regions will permit the application 
of mesh loadings at least 20$ in excess of the 1615 N (363 lb) 
design loading without tip interference — including an allowance 
for modest tooth spacing errors . 


64 



Table 1. Tabulation of Single Tooth Loadings and Deflections 
as a Function of Roll Angle. Loadings were Obtained 
from Fig. 12. [From Ref. 18.] 


Roll 

angle 

rad 

deg 

0.2208 

12.65 

.2271 

13.01 

.2361 

13.53 

.2474 

14.17 

.2609 

14.95 

.2760 

15.81 

.2922 

16.74 

.3084 

17.67 

.3238 

18.55 

.3386 

19.40 

.3529 

20.22 

.3640 

20.86 

.3747 

21.47 

.3876 

22.21 

.4000 

22.92 

.4121 

23.61 

.4237 

24.28 

.4350 

24.92 

.4458 

25.54 

.4557 

26.11 

.4645 

26.61 

.4723 

27.06 

.4793 

27.46 

.4855 

27.82 

.4909 

28.13 

.4956 

28.40 

.4997 

28.63 

.5031 

28.83 

.5059 

28.99 

.5080 

29.11 


Distance from pitch 
point along line of 
action 


Deflection 


20.92 

53.59 

100.3 

158.9 

228.9 

307.2 

363.0 

363.0 

363.0 

363.0 

363.0 

363.0 

363.0 

363.0 

363.0 

363.0 

363.0 

363.0 

339.0 
287.7 

272.2 

201.6 

165.3 

133.2 

105.1 
80.77 
59.57 
41.87 
27.35 
16.46 


5.7412 
1.5026 
2.8184 
4.5072 
6.5641 
8.9797 
1.0875 
1 1200 
1.1554 
1.1930 
1.2345 
1.2600 
1.3035 
1.3501 
1.4004 
1.4537 
1.5090 
1.5680 
1.5132 
1.3438 
1.1713 
1.0058 
8.5060 
7.0455 
5.6937 
4.4691 
3.3726 
2.4210 
1.6135 
9.6342 


2.2603 

E-5 

5.9158 E-5 

1.1096 

E-4 

1.7745 E-4 

2.5843 

E-4 

3.5353 E-4 

4.2813 

E-4 

4.4095 

E-4 

4.5488 

E-4 

4.6970 

E-4 

4.8603 

E-4 

4.9605 

E-4 

5.1320 

E-4 

5.3155 E-4 

5.5135 

E-4 

5.7233 

E-4 

5.9408 

E-4 

6.1733 

E-4 

5.9573 

E-4 

5.2905 

E-4 

4.6113 

E-4 

3.9598 E-4 

3.3488 

E-4 

2.7738 

E-4 

2.2416 

E-4 

1.7595 

E-4 

1.3278 

E-4 

9.5315 E-5 

6.3525 

E-5 

3.7930 

£-5 
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Table 2. Tabulation of Single Tooth Loadings and Profile 

Relief as a Function of Roll Angle for Approximate 
Minimization of Vibration Excitation. 


Roll 

angle 

Distance 

from pitch 

Load 

Profile 

relief 

rad 

deg 

Mu mt c long i me ut 

action 

N 

lb 

mm 

in. 



mm 

in. 





0.2208 

12.65 

-5.979 

-0.2354 

93.05 

20.92 

1.5106 E-2 

5.9473 E-4 

.2271 

13.01 

-5.718 

-.2251 

238.4 

53.59 

1.4178 E-2 

5.5817 E-4 

.2361 

13.53 

-5.342 

-.2103 

446.1 

100.3 

1.2862 E-2 

5.0637 E-4 

.2474 

14.17 

-4.867 

-.1916 

706.8 

158.9 

1.1173 E-2 

4.3988 E-4 

.2609 

14.95 

-4.308 

-.1696 

1018 

228.9 

9.1161 E-3 

3.5890 E-4 

.2760 

15.81 

-3.673 

-.1446 

1366 

307.2 

6.7005 E-3 

2.6380 E-4 

.2922 

16.74 

-2.997 

-.1180 

1615 

363.0 

4.8057 E-3 

1.8920 E-4 

.3084 

17.67 

-2.322 

-.0914 

1615 

363.0 

4.4801 E-3 

1.7638 E-4 

.3238 

18.55 

-1.676 

-.0660 

1615 

363.0 

4.1262 E-3 

1.6245 E-4 

.3386 

19.40 

-1.057 

-.0416 

1615 

363.0 

3.7498 E-3 

1.4763 E-4 

.3529 

20.22 

-.462 

-.0182 

1615 

363.0 

3.3350 E-3 

1.3130 E-4 

.3640 

20.86 

0 

0 

1615 

363.0 

3.0805 E-3 

1.2128 E-4 

.3747 

21.47 

.447 

.0176 

1615 

363.0 

2.6449 E-3 

1.0413 E-4 

.3876 

22.21 

.986 

.0388 

1615 

363.0 

2.1788 E-3 

.8578 E-4 

.4000 

22.92 

1.506 

.0593 

1615 

363.0 

1.6759 E-3 

.6598 E-4 

.4121 

23.61 

2.012 

.0792 

1615 

363.0 

1.1430 E-3 

.4500 E-4 

.4237 

.4350 

24.28 

24.92 

2.497 

2.969 

.0983 

.1169 

1615 

1615 

363.0 

363.0 

.5906 E-3 

0 

.2325 E-4 

0 

.4458 

25.54 

3.419 

.1346 

1508 

339.0 

.5486 E-3 

.2160 E-4 

.4557 

26.11 

3.830 

.1508 

1280 

287.7 

2.2423 E-3 

.8828 E-4 

.4645 

26.61 

4.199 

.1653 

1077 

272.2 

3.9675 E-3 

1.5620 E-4 

.4723 

27.06 

4.526 

.1782 

896.7 

201.6 

5.6223 E-3 

2.2135 E-4 

.4793 

27.46 

4.818 

.1897 

735.3 

165.3 

7.1742 E-3 

2.8245 E-4 

.4855 

27.82 

5.077 

.1999 

592.5 

133.2 

8.6347 E-3 

3.3995 E-4 

.4909 

28.13 

5.304 

.2088 

467.5 

105.1 

9.9865 E-3 

3.9317 E-4 

.4956 

28.40 

5.502 

.2166 

359.3 

80.77 

1.1211 E-2 

4.4138 E-4 

.4997 

28.63 

5.669 

.2232 

265.0 

59.57 

1.2308 E-2 

4.8455 E-4 

.5031 

28.83 

5.812 

.2288 

186.2 

41.87 

1.3259 E-2 

5.2201 E-4 

.5059 

28.99 

5.928 

.2334 

121.7 

27.35 

1.4067 E-2 

5.5380 E-4 

.5080 

29.11 

6.017 

.2369 

73.21 

16.46 

1.4717 E-2 

5.7940 E-4 
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FIG. 13. PROFILE MODIFICATION FOR APPROXIMATE 
MINIMIZATION OF VIBRATION EXCITATION. 
[FROM TOOTH DEFLECTION CURVE GIVEN 
IN REF. 18.] 
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Transmission Error Dynamic Increment Functions 

In the section devoted to the response, we found that the 
dynamic contributions of the static transmission error to the 
mesh and bearing support force spectrum, and to the force on a 
single tooth, all are governed by the functions QjjtD and 

Q j jj(w), where, in general, Q^. (w) represents the dynamic con- 
tribution to the transmission error complex harmonic amplitude 
at mesh i from a unit-amplitude harmonic component of the static 
transmission error at mesh j. Here, either subscript i or j can 
represent either mesh I or mesh II, where in the foregoing, we 
have been concerned mainly with the case i = I. See, for example, 
Eqs . (3.7), (3.8), and (3.11), and Eqs . (3- 3D, (3-33), and 

(3.34). The "transfer functions" Q. . (to) are properties only of 

a J 

the gearing structure — we have designated them "transmission 
error dynamic increment functions" since each function Q. .(to) 
tends to zero as w -*■ 0 . 

For the four-square test apparatus under study, expressions 
for Qjj(D in terms of basic system parameters are given by 

Eqs. (2.51), (2.55), (2.57), and ( 2 . 58 ), and expressions for 
Q x jqCD are given by Eqs. (2.52), ( 2 . 56 ), (2.57), and (2.59). 
Asymptotic low-frequency approximations for Qjj(<D and 

are given by Eqs. (2.62) and (2.65) respectively, and asymptotic 
high-frequency approximations for Q^(u)) and ^^(w) are given 

by Eqs. (2.64) and (2.66) respectively. Functions Qjj(D and 
Ql pj( w ) are, in general, complex. 

The basic system parameters required to evaluate the above 
expressions have been carefully evaluated and are: 

I = .004112 N-m-sec 2 (.0364 lb-in-sec 2 ) 

1 

I = .0008010 N-m-sec 2 (.00709 lb-in-sec 2 ) 

2 

I = .002994 N-m-sec 2 (.0265 lb-in-sec 2 ) 

3 

R = R b = .0417693 m (1.64446 in). 

C = 1.26312 x 10“ 8 m/N (2.21195 x 10” 6 in/lb) 

C JX = 2.10520 x 10“ 9 m/N (3.68659 x 10- 7 in/lb). (5-5) 

Values of the damping constants c_, c „ , c , and c have been 

1 z d a 

chosen by assuming* that the damping coefficients of all four 


*Private communication with Dr. John J. Coy. 
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bearings in the test rig are the same, where the mesh damping 
is included in the bearing damping coefficients. The bearing 
damping is then specified by the fraction of critical damping 
C in the test gear-shaft-slave gear degree-of-freedom on the 

side of the unit opposite to the loading-vane; the loading-vane 
damping coefficient c is specified by assuming a value for the 

fraction of critical damping L, in the slave gear degree-of- 

3 

freedom on the loading-vane side of the unit, which includes the 
slave gear bearing damping and the damping associated with rela- 
tive motion in the loading-vane. Expressions for c and £ are 

13 

readily determined from Eqs . (2.19) and (2.21) as 


r — 

c c + c 

1 j. _ 3 a 

1 

2/k I 3 2/k I 

11 33 

where 


k = 
1 

hH^)- 1 + (c II )" 1 ] , k 3 

from which 

we obtain 


(5 • 6a,b ) 


(5 • 7a,b ) 


c =25 /k I 
1 111 


c = (2C y/k~T ~ ) - c , 
a 333 3 


and 



(5 . 8a,b ) 


(5.9) 


This last expression is a consequence of the assumption of equal 
bearing damping coefficients and the fact that c^ represents the 

damping of two bearings, whereas c and c represent the damping 

2 3 

of only a single bearing each. Assuming* values of = 0.1 and 

C 3 = 0.5 yields, from Eqs. (5.5) through (5-9)> the following 
damping constants: 


c^ = 12.6114 N-m-sec (111.626 Ib-in-sec) 

c = c = 6.30571 N-m-sec (55-8130 lb-in-sec) 

2 3 

c o = 43.5062 N-m-sec ( 385 . 081 Ib-in-sec). (5-10) 

3 . 


*Private communication with Dr. John J. Coy. 
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The above cited expressions for the two transmission error 
dynamic increment functions and Q-j. jj(^) and their low- 

and high-frequency asymptotic approximations have been evaluated 
using the above cited numerical constants, and are plotted in 
Pigs. 14 through 17. Figure 14 displays the magnitude of Q^Cm) 

and the magnitude of the low- and high-frequency asymptotic ap- 
proximations of Qjj-(to), Pig. 15 displays these same metrics for 

the magnitude of Q^. ^(w). Fig. 16 displays the real and imagi- 
nary parts of Qjj( w ) and its asymptotic approximations, and 
Pig. 17 displays the comparable quantities for ^(w). 

The test apparatus shown in Pig. 1 has three degrees of 
freedom; however, one of these must be considered to be degenerate 
due to the absence of an elastic restraint on steady rotation of 
the system. Consequently, we would expdct to see two resonant 
frequencies for such a system. Figures 14 and 15 show one reso- 
nant frequency in the neighborhood of about 3500 Hz, but no 
second resonance is apparent in these figures. This second reso- 
nance does not show up because of the relatively large value 
(C = 0 . 5 ) of damping assumed for relative motion in the loading- 

3 

vane. Choice of a smaller value for £ shows a second resonance 

3 

in the approximate neighborhood of 2000 Hz. 

The nominal running speed of the test apparatus is 10,000 
rpm. Thus, the tooth meshing frequency of the 28 tooth test gears 
is 

f I = ^ q ° ° x 28 = 4666.67 Hz, ( 5 . 11) 

and the tooth meshing frequency of the 35 tooth slave gears is 

f II = ~ 6 o° ~ x 35 = 5833.33 Hz . (5.12) 

Examination of Pig. 14 shows that f T is beyond both system resonant 
frequencies, and is in the region where the high-frequency asymp- 
totic formula for Q^Cw) given by Eq. (2.64) provides quite 

accurate results. Examination of Pig. 15 shows that fjj also is 

beyond both resonant frequencies, but the high-frequency asymp- 
totic formula, Eq. (2.66), provides only marginal accuracy for 
Ql jj(to) in the neighborhood of f^ = 5833 Hz. 

Further examination of Fig. 14 shows that for fj- = 4667 Hz 
and all of its multiples, | (co) | is approximately unity indicat- 
ing that dynamic effects are playing an important part at all of 
the test gear tooth-meshing harmonic frequencies. On the other 
hand, examination of Fig. 15 shows that for f^ = 5833 Hz and 



DYNAMIC INCREMENT FUNCTION 


Q i < W-H n f ;> 




DYNAMIC INCREMENT FUNCTION 
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beyond, | II (co)| is not very different from zero, indicating 

that at the 10,000 rpm testing speed, the tooth-meshing har- 
monics of the static transmission error of the slave gears have 
little effect on the dynamic transmission error of the test 
gears. However, if the test apparatus rpm were reduced to 
about 6000 rpm. Pig. 15 shows that a very substantial magnifi- 
cation of the slave gear static transmission error tooth meshing 
fundamental would be added to the test gear dynamic transmission 
error. See Eq. (3.17). This magnification factor , which is 
about 5.6j provides dramatic illustration of the very substantial 
errors that can be introduced by the common practice of modeling 
gearing systems with more than one mesh by simple uncoupled 
single mesh systems. 

Mesh and Bearing-Support Force Spectrum 

In the remaining portion of this report, we illustrate the 
theory developed in the preceding portions with computations 
made for two sets of test gears: one set possesses perfect in- 

volute profiles except for linear tip relief that starts at the 
highest point of single tooth contact and extends to the tip of 
the tooth; the second set possesses the approximately optimum 
profile modification defined by Eq . (5.4), where the tooth loading 

function W.(x) used is illustrated in Fig. 12. Thus, the teeth of 
J 

this second gear set have additions of material to the perfect 
involute profiles, where the added depths are equal to the values 
of the deflections given in Table 1. Linear interpolation was 
used to define the values of these profile modifications between 
the values given in Table 1 — except for the points on either 
side of the slope discont inuities shown in Fig. 13. Additional 
points were added at the locations of these slope discontinuities 
that were computed by linear extrapolation of pairs of points 
on either side of the discontinuities shown in Fig. 13. These 
two additional points are: 


Roll 

angle 

"Deflection" 

rad 

deg 

mm 

in 

.2872 

16.46 

1.0772 E-2 

4.2409 E-4 

.4408 

25 . 26 

1.5984 E-2 

6.2929 E-4. 


The above values can be thought of as added entries to Table 1. 

The amount of tip relief used in the linear tip relief test 
gear calculations was .01524 mm (.0006 in) which was specified 
at a radial position of the addendum radius minus the tip rounding 
radius of .381 mm (.015 in). 


75 



The first step in determining the mesh and bearing-support 
force spectrum, Eq. (3.11), is computation of the sets of Fourier 
series coefficients a mnJ; and a mnII , n = ±1,±2,**- of the static 

transmission errors of the test gears (mesh I) and those of the 
slave gears (mesh II) for use In Eqs. (3*7) and (3.8). The tooth 
profiles of the slave gears are unmodified involute profiles. The 
sets of static transmission error Fourier series coefficients, 

Eq. (1.5), were computed using the general methods outlined in 
Sec. 1 of this report. Each set of Fourier series coefficients was 

then multiplied by (C-,) -1 times the appropriate dynamic increment 
function as indicated x by Eqs. (3-7) and (3-8) to yield the test 
gear mesh and bearing-support force Fourier series coefficients 
arising from the test gear and slave gear static transmission error 
excitations. The dynamic increment functions used were those 
illustrated in Figs. 14 through 17. The test gear mesh and bearing- 
support force Fourier series coefficients a yjj( n ) and jj(n) 

that arise from the static transmission errors of the test and 
slave gear meshes respectively were then combined as in Eq. (3-11) 
with the d-c torque and base circle radius R of the test gears 

to yield the test gear mesh and bearing-support force (line) spec- 
trum. Frequencies f^ and f^ are the tooth meshing frequencies 
of the test and slave gears respectively. 

Figure 18 shows the (one-sided) line spectrum of the dynamic 
loads experienced by the test gear mesh and bearings that arises 
from the test gear static transmission error, and that is ex- 
pressed by the first line of Eq. (3.11). The spectrum shown in 
Fig. 18 is for the above-mentioned test gear set with linear tip 
relief. Figure 19 shows the comparable spectrum for the test gear 
set with the approximately optimum profile modification shown in 
Fig. 13. The line spacing in each of Figs. 18 and 19 is the 
frequency spacing of the test gear tooth meshing harmonics, which 
is 4666.67 Hz according to Eq. (5.11). Each plot shows a d-c 
component on the ordinate with strength marked by an x and 64 
additional lines, which was the total number of harmonics computed 
by the computer program. The abscissa is measured in Hertz and 
the (logarithmic) ordinate is measured in (lb) 2 . 

The value of the d-c component indicated by the x is the same 
in each spectrum shown in Figs. 18 and 19. This value is 
(x /R) 2 = ( 363) 2 = 131,769 (lb) 2 . However, careful inspection of 

cl 

the remaining lines shows that those of Fig. 18 are typically 
about 10 3 stronger than those of Fig. 19 corresponding to (linear) 

forces in Fig. 19 of about /10“ 3 or 3% to of those in Fig. 18. 

In fact, since the tooth optimum profile modification is supposed 
to exactly compensate for the static transmission error, the 
strengths of all lines in Fig. 19 except the d-c component should 
theoretically be zero. This did not happen exactly because 
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FIG. 18. PREDICTED CONTRIBUTION FROM TEST GEAR MESH STATIC 
TRANSMISSION ERROR TO (ONE-SIDED) FORCE POWER 
SPECTRUM OF TEST GEAR MESH AND BEARING SUPPORTS 
FOR CASE OF LINEAR TIP RELIEF OF TEST GEARS. 
[FIRST LINE OF EQ. (3.11).] 
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FIG. 19. PREDICTED CONTRIBUTION FROM TEST GEAR MESH STATIC 
TRANSMISSION ERROR TO (ONE-SIDED) FORCE POWER 
SPECTRUM OF TEST GEAR MESH AND BEARING SUPPORTS FOR 
CASE OF APPROXIMATELY OPTIMUM PROFILE MODIFICATION 
OF TEST GEARS. [FIRST LINE OF EQ. (3.11).] 
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we modeled the single tooth compliance by a least-squares aubia 
polynomial fit to the compliance values computed from the data 
shown in Table 1, whereas we took the profile modifications 
directly from the numerical values shown in Table 2 — giving rise to 
a small discrepancy between the compliance values and deflection 
values caused by the small "errors" in the cubic polynomial fit. 
However, such force spectrum reductions shown by comparison of 
Figs. 19 and 18 might be typical of what we could hope to achieve 
with the optimum design procedure indicated by Eqs . (4.6a) and 
(4.6b). We would not expect to do this well in practice with the 
approximate optimum design given by Eq. (5.4), which neglects de- 
flections of all teeth except the loaded tooth. 

Figures 20 and 21 are plots of the contributions of the second 
line in Eq . (3.11) to the mesh and bearing support force spectra. 
Figure 20 provides that contribution in the case of the test gears 
with linear tip relief, and Fig. 21 provides the comparable con- 
tribution in the case of the test gears with the approximately 
optimum profile modifications prescribed by Table 2.- The line spac- 
ing in each of Figs. 20 and 21 is the frequency spacing of the slave 
gear tooth meshing harmonics which is 5833-33 Hz, according to 
Eq . (5.12). The abscissa of Figs. 20 and 21 is measured in units 

of Hertz and the (logarithmic) ordinate is measured in units of 
(lb) 2 . 

The line spectra shown in Figs. 20 and 21 possess no d-c com- 
ponent. This lack of a d-c component is a result of the fact that 
Qj jj(oj) = 0 at w = 0 as can be seen from Fig. 15 — which results 

in jj(0) = 0 according to Eq. (3.8). 

According to Eqs. (3.8) and (3.11), the strengths of the lines 
in Figs. 20 and 21 should depend on the time-average value of the 

compliance of mesh I, C ^ , the transmission error "cross" dynamic 

increment function shown in Fig. 15, and the Fourier series 

coefficients a TX of the static transmission error of (slave) gear 
mnii 

mesh II. Since trie slave gears had no profile modifications, 
their static transmission error arises only from their elastic de- 
formations. The Fourier series coefficients of the elastic de- 
formation component of the static transmission error is given by 
Eq. (1.8) — they are determined essentially by the Fourier series 
coefficients a (p/K) n > n = 0,±1,±2,’** of the instantaneous recipro- 
cal of the total stiffness of mesh II. For the slave gears of 
mesh II, these Fourier series coefficients were taken to be one- 
sixth of those for the test gears of mesh I (because the width of 
the slave gears is six times that of the test gears). Since the 
roll distance L of the test gears with linear tip relief differs 
somewhat from the design roll distance L for the test gears with 
(approximately) optimum profile modifications, the temporal behavior 
of the test gear instantaneous total mesh stiffness differs slightly 
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FIG. 20. PREDICTED CONTRIBUTION FROM SLAVE GEAR MESH STATIC 
TRANSMISSION ERROR TO (ONE-SIDED) FORCE POWER 
SPECTRUM OF TEST GEAR MESH AND BEARING SUPPORTS 
FOR CASE OF LINEAR TIP RELIEF OF TEST GEARS. 
[SECOND LINE OF EQ. (3.11).] 
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FIG. 21. PREDICTED CONTRIBUTION FROM SLAVE GEAR MESH STATIC 
TRANSMISSION ERROR TO (ONE-SIDED) FORCE POWER 
SPECTRUM OF TEST GEAR MESH AND BEARING SUPPORTS FOR 
CASE OF APPROXIMATELY OPTIMUM PROFILE MODIFICATION 
OF TEST GEARS. [SECOND LINE OF EQ. (3.11).] 



for these two cases — hence, our predictions of the slave gear 
contributions to the test gear mesh and bearing-support force 
spectra shown in Figs. 20 and 21 also differ slightly for these 
two cases. 


From Eqs . (3-7), (3*8), and (3.11), we see that the trans- 
mission error "auto" dynamic increment function Qjj(w) of Fig. 14 

is the structural transfer function that governs the force spectrum 
component given by the first line in Eq. (3.11), whereas the trans- 
mission error "cross" dynamic increment function ^(w) of Pig. 15 

is the structural transfer function that governs the force spectrum 
component given by the second line in Eq. (3.11). In the frequency 
range beyond 4666 Hz where all test gear tooth-meshing harmonics 
fall, we see from Fig. 14 that [ (oj ) [ = 1, whereas in the fre- 
quency range beyond 5833 Hz where all slave gear tooth meshing 
harmonics fall, we see from Fig. 15 that j jj(w)| << 1. This 

difference in the structural transfer functions is the main reason 
why the slave gear contribution of Fig. 20 to the test gear force 
spectrum is much less than the test gear "auto" contribution shown 
in Fig. 18 for the linear tip relief case. However, in the case 
of the approximately optimum profile modification, the first two 
slave gear tooth meshing harmonic contributions shown in Fig. 21 
are stronger than those from the test gear shown in Fig. 19 be- 
cause of the reduction of the test gear contributions arising from 
the approximately optimum profile modification. 


Dynamic Force on a Single Tooth 

The dynamic increment 6W^j(t) of the force history on a single 
tooth given by Eq. (3-31) can be decomposed into two components 


«w+ 3 (t) - aw+. l( t) + M+jii (t ) . 

(5.13) 

where we have defined 


«w+ 31 (t) 4 K+jCt) «4 IZ (t) 

(5.14) 

and 


= K Sj (t) 6 c ml II (t) J 

(5.15) 


where the notation on the right-hand sides of Eqs. (5.14) and (5.15) 

is the same as that on the right-hand side of Eq. (3.31). That is, 

+ 

K . (t) is the instantaneous stiffness of a typical pair of test 

^ J 4- 4- 

gear teeth, and S£ m jj(t) and 65 ^ jj(t) are the dynamic increments 

of the transmission error of the test gear mesh arising from the 
test gears and slave gears as given by Eqs. (3-33) and (3*34) re- 
spectively. B'rom Eqs. (3.29), (3-31), and (5.13) through (5.15), 
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we can represent the dynamic force W--.(t) on a single tooth of 
the test gear mesh as J 


W 


Ij 


(t) = [w+ j 0 (t) + 




>] 




(5.16) 


The first two terms within the brackets in the right-hand side of 
Eq. (5.16) constitute tht total contribution to the single tooth 

force W T ,(t) arising from' the static transmission error excitation 


IJ 


rt 


by the test gear mesh, whereas the third term 6W.^. .j. ^ ( t ) in the 
right-hand side of Eq. (5.16) is the contribution to Wj. (t) from 
the static transmission error excitation by the slave gear mesh. 


Figure 22 shows the predicted total force contribution 

[Wj. j 0 ^ ^ +( ^Ij I ^ ^ on a t yp lcal to °6h of the driving test gear 

arising from the static transmission error excitation by the test 
gear mesh. The rotational speed of the NASA Lewis four-square 
test apparatus used in the computation was 10,000 rpm. The test 
gear mesh carries a mean loading x /R of 1615 N (363 lbs). The 

OL 

force history shown in Fig. 22 is for the test gears with linear 
tip relief of .01524 mm (.0006 in) as described earlier in this 
report . 


The ordinate of Fig. 22 is force in units of pounds, and the 
abscissa is time in seconds. The time origin t = 0 is located at 
the pitch point which corresponds to the center position of the 
time axis. Contact initiation occurs at the base of the tooth at 

t = - 1 . 57 x 10“ 4 sec and terminates at the tip at t = 1.57 x 10 -4 sec. 

At the instant of contact initiation, the tooth abruptly picks 
up a loading of 98.5 lb as is indicated in Fig. 22. This abrupt 
loading of the base occurs because the tip of the tooth of the 
mating gear has too little modification to allow a smooth onset of 
loading of the tooth base at contact initiation. The loading on 
the tooth thereupon increases in an approximately linear fashion 
from the instant of contact initiation through the interval of 
two-tooth contact until the instant where the discontinuity in 
slope occurs and where the rate of loading increase abruptly di- 
minishes. Beyond this instant, where the force on the tooth is 
approximately a constant value of 338 lb, is the region of single 
tooth pair contact. The loading on the tooth thereupon decreases 
in an approximately linear fashion again until contact terminates 
at the tip at t = 1.57 x 10 -4 sec. At this instant, the loading on 
the tooth drops abruptly from 98.5 lb to zero as contact is lost. 

Figure 22 indicates that the maximum loading carried by the 
tooth is 338 lb. This value is 25 lb less than the static loading 
of 363 lb carried by the tooth within the region of single tooth 
pair contact under static loading conditions. The reason for this 
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FIG. 22. PREDICTED FORCE HISTORY ON A TYPICAL DRIVING TEST GEAR 
TOOTH WITH LINEAR TIP RELIEF. ROTATIONAL SPEED IS 
10,000 RPM. MESH CARRIES A NOMINAL LOADING OF 
1615N (363 LB). TIP RELIEF IS .01524 MM (.0006 IN). 
FIGURE EXCLUDES EFFECTS OF DYNAMIC EXCITATION FROM 
SLAVE GEARS. 
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■ K 


25 lb maximum loading decrease is that at the instant of maximum 
loading of 338 lb shown in Fig. 22, the dynamic loading increment 

Swt. T (t) of Eq. (5.14) is -25 lb. This loading decrease arises 

from the dynamic transmission error increment | 5? m jj(b) in Eq . 

(5.1^) (which is negative in this case). 

Physically, this decrease of 25 lb of maximum tooth loading 
arises from the undermodification of the tooth tips that results 
in the abrupt loading increments of 98.5 lb at contact initiation 
and termination as shown in Fig. 22. This "extra" loading carried 
in the two tooth-pair contact regions results in a reduced tooth 
loading within the region of single tooth-pair contact, which occurs 
because at a test speed of 10,000 rpm the inertia associated with 
the test gear masses prevents the teeth from picking up their full 
static loading of 363 lb within the single tooth-pair contact 
interval . 

*f* *4* 

Figure 23 shows the force history [W^ q ( t )+6W^j ^ ( t ) ] on a 

typical tooth of the driving test gear that is comparable to Fig. 

22 — except that the force history shown in Fig. 23 was computed 
for test gears with the profile relief values shown in Table 2 
(which was displayed earlier in this section) . Force values 

-I* *4« 

[W-j-j q (t ) + 6W.j. . i(t ) ] shown in Fig. 23 arise from the static trans- 
mission error excitation by the test gear mesh only, as can be 
seen from Eqs. (3.29), (3.30), (3.33), (5.13), and (5.14). For 
the case shown in Fig. 23, the rotational speed of the four-square 
test apparatus was 10,000 rpm and the mean loading x /R carried by 

the test gears was 1615 N (363 lb), which are the same values as 
for the case shown in Fig. 22. 

Contact initiation occurs at the base of the tooth at a value 

of t = -1.41 x 10 -4 sec instead of t = -1.57 x 10 -4 sec as in 
Fig. 22. This reduction in |t| at contact initiation (and contact 
termination) shown in Fig. 23 occurs because of the slightly 
smaller design value of roll distance of L = 1.22971 cm (.48414 in) 
of the test gears of Fig. 23 [in comparison with 

L = 1.37450 cm (.5411406 in) for the test gears of Fig. 22]. In 
contrast to the case shown in Fig. 22, the tooth loading at the 
instant of contact initiation of t = -1.4l x 10 -4 sec and contact 

termination t = 1.4l x 10 -4 sec is zero for the case shown in 
Fig. 23 . The initial value of 17.6 lb shown in Fig. 23 occurs at 

a value of |t| somewhat smaller than |t) = 1.41 x 10” 4 sec. The 

loading history between |t| = 1.4l x 10 -4 sec and the first value 
of 17.6 lb shown in Fig. 23 is a linear extrapolation of the his- 
tory shown in Fig. 23. This small region is missing from the fig- 
ure because of an artifact of the computer plot routine used to 
plot Fig. 23 . 
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PREDICTED FORCE HISTORY ON A TYPICAL DRIVING TEST 
GEAR TOOTH WITH "OPTIMUM" PROFILE RELIEF. ROTATIONAL 
SPEED IS 10,000 RPM. MESH CARRIES A NOMINAL LOADING 
OF 1615N (363 LB). PROFILE RELIEF VALUES ARE GIVEN 
IN TABLE 2. FIGURE EXCLUDES EFFECTS OF DYNAMIC 
EXCITATION FROM SLAVE GEARS. 
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For the case shown in Fig. 23, the single tooth loading pro- 
gresses linearly from zero at t = -1.41 x io -4 sec to a value of 
about 363 lb at the instant where the region of single tooth pair 
contact begins. The loading during the region of single tooth 
pair contact is almost a constant value growing from about 363 lb 
to 366 lb and back to 363 lb at the end of the region of single 
tooth pair contact. During the subsequent region of two tooth pair 
contact, the loading decreases linearly from 363 lb to zero where 
contact terminates at the tip. The fluctuation of about 3 lb 
within the region of single tooth pair contact occurs because the 
cubic polynomial representation of single tooth compliance that we 
have used in the computation does not agree exactly with the stiff- 
ness values computed from the deflections shown in Table 1 — which 
were used to compute the profile relief values shown in Table 2 
(as noted earlier) . 


Exact implementation of the optimum design procedure results 
in zero static transmission error ? m j(b) for the optimally de- 
signed teeth. Thus, for such idealized cases, all Fourier series 
coefficients a T of the static transmission error are zero: hence, 

the dynamic increment ^ (t ) also is identically zero for these 
cases, as we may see from Eq . (3-33). It then follows from Eq. 

4 * 

(5.1^) that the dynamic loading increment <5W^ . ^(t) also is iden- 
tically zero in these idealized cases so that, for a perfectly de- 


signed profile modification. Fig. 23 should be a plot of W-^. Q (t) 

which is the loading history experienced by a typical tooth under 
quasi-static running conditions when inertial effects are negligi- 
ble. This loading history is the prescribed single-tooth loading 
function used in the optimal design process — which, for the pres- 
ent case, is shown in Fig. 12. Thus, except for the scaling of 
the independent variable from roll distance to time, the loading 
functions in Figs. 12 and 23 should be identical. Examination of 
the two figures shows that they are very nearly identical except 
for the approximately 3 lb fluctuation in Fig. 23 within the 
region of single tooth pair contact that arises from the "error" 
in the cubic polynomial representation of the single tooth com- 
pliance data computed from Table 1. 


Comparison of the tooth loading histories shown in Figs. 22 
and 23 shows that the maximum value of single tooth loading for the 
tooth designed for minimum vibration excitation is some 28 lbs 
larger than the maximum single tooth loading for the teeth de- 
signed with linear tip relief (with somewhat undercorrected tips). 
The reason for this result was implied earlier — namely, the abrupt 
(step) loading at contact initiation (and termination) of 98.5 lb 
in the linear tip relief case (that resulted from "undercorrection" 
of the tips) permits the region of two-tooth pair contact to effec- 
tively carry some of the loading that would normally be carried 
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within the single tooth-pair contact region at lower rotational 
speeds. This reduction in peak' single tooth loading of 28 lbs 
is a consequence of the rotary inertias associated with the test 
gears, and would not be observed at substantially lower rotational 
speeds. Furthermore , excitation of other parts of a gearing system 
arising from such undercorrected profiles oould — and often will — 
result in loading increases of such other parts of the system as 
we shall next show. 


f 

In order to obtain the total force W x ( t ) on a typical test 

i 

gear tooth, the dynamic force increment ^ w ijji(t ) mus t be added 

to the quantity [wt, ( t ) + <5wt . T (t ) ] plotted in Pigs. 22 and 23. 

This fact can be seen from Eq. (5-16), where SW-j-jjj(t) represents 

the dynamic force increment on a typical tooth j in the test gear 
mesh (I) from the static transmission error excitation arising from 

the slave gear mesh (II). The dynamic force increments ( t ) 

computed for the case of linear tip relief, and then for the case 
of the profile relief values of Table 2, are shown in Pigs. 24 and 
25 respectively. These two figures differ for the same reason that 
the spectra shown in Pigs. 20 and 21 differ. That is, in each com- 

4 * 

putation of SW T . TT (t), we require the dynamic transmission error 


. t 


IjH 


increment yj(t) according to Eq. (5.15) , which in turn is de- 

termined from the Fourier series coefficients a TT of the static 

mnll 

transmission error of the slave gears as is indicated by Eq. (3-34). 
These Fourier series coefficients a mn jj are determined essentially 

by the Fourier series coefficients °f the instantaneous 

reciprocal of the total mesh stiffness of the slave gear mesh, 
which were taken to be one-sixth of those for the test gear as 
indicated earlier. However, the roll distance L over which tooth 
contact takes place differs for the linear tip relief case in com- 
parison with the profile relief case tabulated in Table 2; hence, 
the temporal behavior of the instantaneous total mesh stiffness 
is different for the two cases, as are the Fourier series coeffi- 
cients of the instantaneous reciprocal total mesh stiffness. 


Figures 2<4 and 25 both are dominated by a single sinusoidal 

component with a period of about T = 2x8.4xio -5 sec ~ 17*10 -5 sec, 
as is readily apparent from the time-axis markings in Fig. 25- 
The frequency of this sinusoidal component is 

(1/T) « (17 x 10 -5 ) ~ 5880 Hz. Since, according to Eq . (5.12), 

the slave gear tooth meshing fundamental frequency is 5833 Hz, 
it is clear from the above approximate calculation that the domi- 
nant sinusoidal component in Figs. 24 and 25 arises from the tooth 
meshing fundamental frequency of the slave gears. This conclusion 
is consistent with Eqs . (3-3*0 and (5.15), and the fact that the 


88 


FORCE INCREMENT, lb 


7.07ET+00 


4.31 E+00 


1 .55E+00 


-1 .21E+00 


-3.97E+00 


*6 . 72E+00 



i l l 



03 

CjJ 

Od 

CD 


» 

■ 

■ 

• 

■ 

m 

CXI 

4k 



4k 

cn 


CjJ 

A 

-0. 

CD 

*\I 

m 

i 

m 

t 

n 

i 

n 

■ 

m 

t 

n 

t 

• 

es 

i 

<s 

C3 

I 

<3 

C3 

i 

Q 

4* 

CJl 

U1 

Ul 

tn 

4w 


TIME, sec 


FIG. 24. PREDICTED FORCE INCREMENT ON A TYPICAL DRIVING TEST 
GEAR TOOTH WITH LINEAR TIP RELIEF THAT RESULTS FROM 
DYNAMIC EXCITATION PROVIDED BY THE SLAVE GEAR MESH. 
ROTATIONAL SPEED IS 10,000 RPM. TEST GEAR MESH 
CARRIES A NOMINAL LOADING OF 1615N (363 LB). ORIGIN 
OF TIME AXIS IS UNRELATED TO THAT OF FIG. 22. 
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FIG. 25. PREDICTED FORCE INCREMENT ON A TYPICAL DRIVING TEST 
GEAR TOOTH WITH "OPTIMUM" PROFILE RELIEF THAT 
RESULTS FROM DYNAMIC EXCITATION PROVIDED BY THE 
SLAVE GEAR MESH. ROTATIONAL SPEED IS 10,000 RPM. 
TEST GEAR MESH CARRIES A NOMINAL LOADING OF 
1615N (363 LB). ORIGIN OF TIME AXIS IS UNRELATED 
TO THAT OF FIG. 23. 
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transmission error "cross" dynamic increment function jj( w ) 

is very nearly zero for all multiples of the slave gear tooth 
meshing fundamental frequency except for the first which occurs 
at 5833 Hz as we can see from Pig. 15 . 

The additional loading of about 5 to 7 lb added to the test 
gear teeth as a result of the slave gear mesh static transmission 
error excitation is clearly insignificant. However, this situa- 
tion is not necessarily typical. For example, if the test appara- 
tus rotational speed of 10,000 rpm were reduced to 6,000 rpm, 
then the slave gear tooth meshing fundamental frequency would be 
reduced to 0.6 x 5833.33 = 3,500 Hz according to Eq. (5.12). From 
Fig. 15, we see that for this speed reduction, | Q-j- would be 

increased from about |Qj jjl ~ 0.43, which is its approximate 

value at 5833 Hz to I Qj l ~ 5.6 which is its approximate value 

at 3500 Hz. The fractional increase in | Q_ is (5-6/0.43) z 13- 

Consequent ly , at 6,000 rpm running speed, we could expect a maxi- 
mum loading increment from the slave gear mesh of about 
13 x 6 = 78 lb added to the test gear tooth loading from the dy- 
namic excitation provided by the slave gears. This value of 78 lb 
is no longer an insignificant fraction of the 366 to 338 lb maxi- 
mum loading seen in Figs. 22 and 23 that arises from the test gear 
dynamics alone. The major conclusion to be drawn from this simple 
exercise is that a significant loading increment on the teeth of 
one mesh can occur as a result of the dynamic excitation arising 
from another coupled mesh — in this case the slave gear mesh at 
6 > 000 rpm running speed. To avoid the potential for such effects , 
it is desirable to minimize the static transmission error excita- 
tion from every mesh in a system. Such minimization is accomplished 
using the profile modification design equations (4.6a) and (4.6b). 

Since there is no fixed phase relation between the teeth on 
the slave gears and those on the test gears, the time origin in 
Figs. 24 and 25 must be regarded as completely arbitrary relative 
to that in Figs. 22 and 23. Hence, the loading shown in Fig. 24 
has not been superimposed on that shown in Fig. 22, and the load- 
ing shown in Fig. 25 has not been superimposed on that shown in 
Fig. 23. 


Root Stresses on a Single Tooth 

From the loading histories shown in Figs. 22 through 25, tooth 
root stresses can be calculated using Eq . (3.37) and the material 
derived in Appendix B of this report. Such computations were car- 
ried out for the test gears in Fig. 1 using the following parameters 

test gear tooth thickness, t* = 0.635 cm (0.25 in) (5.17) 

test gear tooth root fillet radius, r = 0.1016 cm (0.04 in) ( 5 . 18 ) 
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Heywood parameter, v = 0.25 (5-19) 
base radius included tooth angle, B = 3.98222 deg (5.20) 
A = 4.0909406 cm (I .6106065 in) (5-21) 
h = 0.56958151 cm (0.22424469 in), (5.22) 

0 


where the last four parameters are defined in Appendix B. 

The root stress histories were computed by Eqs . (3.37) and 

(B.l) (which are the same) using the methodology described in 
Appendix B. The stress histories computed from the loading history 

W.(t) in Pig. 22 are shown in Pigs. 26 and 27, which display the 
J 

root stress histories for the driven and driving gears respectively. 
These two stress histories are almost mirror images of each other — 
their lack of perfect symmetry arises from the lack of symmetry 
of the loading history in Fig. 22 about the pitch point. This 
latter lack of symmetry arises from the dynamic increment of the 
test gear pair response, which is not symmetric about the pitch 
point. The stress histories shown in Figs. 26 and 27 are for test 
gears with the linear tip relief as described earlier — these 
stress histories do not include the stress increments caused by 
the static transmission error excitation from the slave gears. 

The root stress histories computed from the loading history 

W.(t) in Fig. 23 are shown in Figs. 28 and 29, which display the 
J 

stress histories for the driven and driving gears respectively. 

These latter two stress histories are almost perfect mirror images 
of each other — their slight lack of symmetry arises from a very 
small lack of perfect symmetry about the pitch point of the load- 
ing history shown in Fig. 23. The stress histories shown in Figs. 

28 and 29 are for test gears with the "optimum" profile relief 
described earlier. Stress increments caused by the static trans- 
mission error excitation from the slave gears are not included in 
Figs. 28 and 29. Since the tooth loadings (shown in Fig. 22) on 
the test gears with linear tip relief generally are less (for 
these particular running conditions) than the tooth loadings (shown 
in Fig. 23) on the test gears with "optimum" profile relief, we 
see that the root stresses shown in Figs'. 26 and 27 generally run 
less than those shown in Figs. 28 and 29. In all four stress 
histories shown in Figs. 26 through 29, the maximum root stress 
occurs at the instant of the highest point of single tooth con- 
tact (nearest to the tooth tip), as expected. 
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PREDICTED ROOT STRESS HISTORY ON A TYPICAL TOOTH 
OF THE DRIVEN TEST GEAR. LINEAR TIP RELIEF CASE 
STRESSES ARE PREDICTED FROM FORCE HISTORY SHOWN 
IN FIG. 22. FIGURE EXCLUDES STRESS INCREMENT 
ARISING FROM DYNAMIC EXCITATION BY SLAVE GEARS. 
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FIG. 27. PREDICTED ROOT STRESS HISTORY ON A TYPICAL TOOTH 

OF THE DRIVING TEST GEAR. LINEAR TIP RELIEF CASE. 
STRESSES ARE PREDICTED FROM FORCE HISTORY SHOWN 
IN FIG. 22. FIGURE EXCLUDES STRESS INCREMENT 
ARISING FROM DYNAMIC EXCITATION BY SLAVE GEARS. 
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FIG. 28. PREDICTED ROOT STRESS HISTORY ON A TYPICAL TOOTH 

OF THE DRIVEN TEST GEAR. "OPTIMUM" PROFILE RELIEF 
CASE. STRESSES ARE PREDICTED FROM FORCE HISTORY 
SHOWN IN FIG. 23. FIGURE EXCLUDES STRESS INCREMENT 
ARISING FROM DYNAMIC EXCITATION BY SLAVE GEARS. 


95 


4 1 E-04 




ROOT STRESS, psi 


4.37E+04 


3.53E+04 


2.G8E+04 


1 .83E+04 


9.8SE+03 


1 .39E+03 




00 

ru 

ru 

00 

• 

. 

■ 

• 

• 


-tv 

00 

00 

4k 


b) 



to 

PI 

n 

i 

rn 

■ 

n 

i 

PI 

1 

1 

Q 

1 

Q 

1 

CS) 

i 

Q 

1 

Q 

4k 

in 

U1 

Ul 

in 


TIME , sec 


FIG. 29. PREDICTED ROOT STRESS HISTORY ON A TYPICAL TOOTH 
OF THE DRIVING TEST GEAR. "OPTIMUM" PROFILE 
RELIEF CASE. STRESSES ARE PREDICTED FROM FORCE 
HISTORY SHOWN IN FIG. 23. FIGURE EXCLUDES STRESS 
INCREMENT ARISING FROM DYNAMIC EXCITATION BY 
SLAVE GEARS. 
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APPENDIX A 


DERIVATION OF EXPRESSION FOR EFFECTIVE 

LOCAL TOOTH PAIR STIFFNESS 

Let and denote, respectively, the elastic 

<3 (J 

deformations of the mating teeth j on meshing gears (1) and 
(2). Let 


A U u> ♦ „«> 


(A. 1) 


denote the combined deformation of the two mating teeth, where 

u^ 2 \ and u. are defined as positive when they are 
J 3 J 

"equivalent" to removal of material from the tooth surfaces 
[ 4 J . Let u.j, denote the component of the combined deformation 

arising from the force transmitted by tooth pair j f . Then 


u . 
J 

U J 


can be expressed as the superposition 


u 


i -L 


u . . f 
33 


(A. 2) 


where the summation over j * includes all tooth pairs j ' in con- 
tact at the particular instant under consideration. Let c jj' 

denote the deformation of the contact point of tooth pair j 
arising from a unit force transmitted by tooth pair j ' as de- 
fined earlier,* and let W., denote the force transmitted by 
tooth pair j ' . Then, J 


u . . , = c . . , W . , 

J J ’ 3 J J ’ 


(A. 3) 


hence, substituting Eq. (A. 3) into Eq. (A. 2), we have 


u j “ S °jj ■ w j 1 • 


(A. 4) 


Let 


for] 


be the inverse of the matrix • Then the 


force W. -transmitted by tooth pair j can be expressed as 
J 


W. = k . , f u . f , 
0 33' 3' 

j 


(A. 5) 


*The quantity 

0 = 0 U) + ^2) 


applies to the tooth pair; thus, we have 

( 1 ) „( 2 ) 


- where c^., and c^., are as defined in 

J«J tiJ J J <JJ J J 

Figure 2. 
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which is the inversion of Eq. (A. 4). Let W denote the total 
force transmitted by all pairs of meshing teeth - i.e., 

W = £ W, 

d J 



(A. 6) 


where Eq. (A. 5) has been used in going to the second line and 
the order of summation has been reversed in going to the third 
line. Interchanging the roles of j and j’ in the last line of 
Eq. (A. 6) , we have 


VJ * 

J 


£ 

j ’ 


j ' j 


u . 
J 


“j • 


(A.T) 


where we have defined the effective local tooth pair stiffness 
as 

A 


K S1 * £ k 

3 j ’ J ' J 


(A. 8) 


Equation (A. 8) expresses K q . in terms of the elements k. . , 

r OJ J J 

of the inverse of the matrix c . r . However, the matrix 

' I J J J 


hr] 


does not have to be inverted to compute K^j . Let W. , 


be the loading distribution in Eq. (A.4_) that yields a constant 

value u. = u independent of j for all u. — i.e., 

J J 


u 


- s <\jj' w 


(A. 9) 


Then from Eq . (A. 5 ), we have in this case 
w, = £ kjj, u 


i » 


= U 


£ k . . , 

y 33 


= u S Vj * 


(A. 10) 
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where the last line in Eq. (A. 10) is a consequence of Maxwell’s 
reciprocal theorem [p. 9 of Reference 15]. Comparing Eqs. (A. 8) 
and (A. 10), we see that 




(A. 11) 


that is, the solution W., to the matrix equation (A. 4) that 

J 

yields a constant deformation u^ = u yields the effective local 

tooth pair stiffness K Q . when W. is divided by u. The above 

derivation is the matrix counterpart to the derivation contained 
in Appendix B of Reference 4. 


Finally, we note that the effective tooth pair stiffness 
K Sj . depends on the position of the point of aontact on the faces 

of the teeth of tooth pair j . When this position is expressed as 
a function of the coordinate z, Kgj becomes the (effective) local 

tooth pair stiffness Kg(z) encountered in the first section of 

this report. Coordinate z is defined in terms of the involute 
roll angle e of the tooth under consideration by 

z = R^e sincf) + c , (A. 12) 


where R^ is the base circle radius, cp is the pressure angle, and 

the constant c is chosen so that the origin of the coordinate z 
is placed at the exact midpoint of the range D of z where tooth 
contact takes place. (Thus, except in the case of identical 
meshing gears, the origin of z is not generally located at the 
pitch point . ) 
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APPENDIX B 


GEOMETRY REQUIRED FOR EVALUATION OF ROOT STRESS FORMULA PARAMETERS 

The modified Heywood formula given by Cornell [8] for the root 
fillet tensile stress in a gear tooth is 



where o is the root fillet tensile stress at the location indicated 
in Fig. B.l, t' is the tooth thickness measured parallel to the gear 
axis, W. is the instantaneous force normal to the tooth surface 
J 1 

that is transmitted by the tooth, v * ^ according to Cornell [8] and 

Heywood [9], and the remaining quantities in Eq. (B.l) are defined 
in Fig. B.l. 


Angle y in Fig. B.l defines the point where the root fillet 
s 

tensile stress is calculated by Eq. (B.l). Cornell [8] provides an 

equation for the value of y that locates the position of the maximum 

s 

root stress. Cornell's equation is 
tany 


(1+0.16A® • 7 )A ± 


sl+1 B. (4+0.4l6A?' 7 )-(i+0.0l6A?- 7 )A.tan<j>’ , 
i i 3 i i w ’ 


(B . 2 ) 


where 


A ± = + 2(l-cosy sl ) , 


(B.3) 


and 


B 1 = + sin ^si ■ 


(B.l)) 


where h and £ are defined in Fig. B.l, and subscripts i and i+1 
o o 

denote that the transcendental equation (B.2) can be solved itera- 
tively for y with i counting the steps in the iteration procedure. 

o 

Once the angle y is determined, the dimension h shown in 

s s 

Fig. B.l also is determined by the formula 


100 









9 


h = h + 2r(l-cosy ) 
so s 


as one can readily show from Fig. B.l. The quantities t’, r, v, and 
h in Eqs. (B.l) through (B.4) are independent of the point of 

application of loading W. and the angle y , and therefore may be 

J s 

regarded as constants. Equations (3-23) and (3-29) of this report 
determine the magnitude of the load Wy . = W. transmitted by a 
generic tooth j of mesh I. ^ J 

To utilize Eqs. (B.l) through (B.4) in the present analysis, 
we must find expressions for the remaining quantities £ g , £ , 

h^, and £ in terms of fixed parameters and a variable readily 

related to time t. Each of the quantities <j>^, £ g , £ , h^, and £ 

is directly dependent on the instantaneous position of tooth pair 
contact — i.e., the instantaneous position of the load W. which 

changes with time as the gears rotate. The "temporal" variable 
that we shall use here to identify with time is the roll angle £ 
of tooth j which is illustrated in Fig. B.2. Also shown there is 
the angle 6 between the tooth centerline and the intersection of 
the involute (active) tooth surface with the base circle. Several 
other quantities required in the present analysis also are defined 
in Figs. B.l and B.2. 


We shall now derive the following expressions for the above 
"time-dependent" parameters: 


^W 

e-6 

(B.5) 

h Tr = 
W 

2R b (l+e 2 ) J5 siny 

(B.6) 

£ = 
s 

£'+rsiny +R, [ (1+e 2 ^cosy-cosB] 
D SO 

(B . 7 ) 

£ ' = 

s 

£ -R, (1+e 2 ) 35 sinytan(e-g ) 

S D 

(B . 8 ) 

£ = 
0 

R^ (1+e 2 ) Js [cosy-sinytan(.E-8 ) ]-£ 

(B.9) 

where 



Y = 

B-e + arctan e. 

(B.10) 


When Eq. (B.10) is substituted into the right-hand sides of Eqs. (B.6) 
through (B.9), the right-hand sides of Eqs. (B.5) through (B.9) all 
are dependent on only fixed parameters except for the time- 
dependent roll angle e. 


To show the validity of Eq. 

that 


(B.5), we note first from Fig. B.l 


4> 


w 



6 . 


(B.ll) 
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From Fig. B.2, we see that 

Y = 6 - (e-a) = g - e + a (B.12) 

and, also, that 

a + 5 + ^ = it 
or 

6 = j - a. (B.13) 

Substituting Eqs . (B.12) and (B.13) into Eq. (B.ll) yields Eq. (B.5) 
To show the validity of Eq. (B.6), we see from Fig. B.l that 

h W ’ 2R w slnY - (B.14) 

Prom Pig. B.2, we further see that 

R w - R b + (R b E)2 ; 


hence , 



( 1+e 2 )** . 


(B.15) 


Substituting Eq. (B.15) into Eq. (B.l4) yields Eq. (B.6). 

Turning to Eq. (B.7) s we may see from Fig. B.l that 

= R^cosy - R^cosg + . (B.l6) 

From Fig. B.l, we further see that 

£ b = £ b + rsln Y s • (B.17) 

Substituting Eqs. (B.15) and (B.17) into Eq. (B.l6) yields Eq. (B.7) 

Next, from Fig. B.l, we note that 

tu, 

*s = £ s - "T tan K * (B.18) 

Substituting Eqs. (B.5) and (B.6) into Eq. (B.l8) yields Eq. (B.8). 
To show the validity of Eq. (B.9), we note from Fig. B.l that 
h w 

£ 0 = R W cos Y " — tan <t>w - A m . (B.19) 
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Substituting Eqs . (B.5), (B.6), and (B.15) into Eq. (B.19) yields 
Eq. (B. 9 ) . 

Finally, to show the validity of Eq. (B.10), we note from 
Fig. B.2 that 

R b e 

tana = -= — = e. (B.20) 

n b 

Substituting a = arctan £ into Eq. (B.12) yields Eq. (B.10). 
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